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We explore the ground states and quantum phase transitions of two-dimensional, spin S = 1/2, 
antiferromagnets by generalizing lattice models and duality transforms introduced by Sachdev and 
Jalabert (Mod. Phys. Lett. B 4, 1043 (1990)). The 'minimal' model for square lattice antiferro- 
magnets is a lattice discretization of the quantum non-linear sigma model, along with Berry phases 
which impose quantization of spin. With full SU(2) spin rotation invariance, we find a magnetically 
ordered ground state with Neel order at weak coupling, and a confining paramagnetic ground state 
with bond charge {e.g. spin Peierls) order at strong coupling. We study the mechanisms by which 
these two states are connected in intermediate coupling. We extend the minimal model to study 
different routes to fractionalization and deconfinement in the ground state, and also generalize it to 
cases with a uniaxial anisotropy (the spin symmetry group is then U(l)). For the latter systems, 
fractionalization can appear by the pairing of vortices in the staggered spin order in the easy-plane; 
however, we argue that this route does not survive the restoration of SU(2) spin symmetry. For 
SU(2) invariant systems we study a separate route to fractionalization associated with the Higgs 
phase of a complex boson measuring non-collinear, spiral spin correlations: we present phase di- 
agrams displaying competition between magnetic order, bond charge order, and fractionalization, 
and discuss the nature of the quantum transitions between the various states. A strong check on our 
methods is provided by their application to S — 1/2 frustrated antiferromagnets in one dimension: 
here, our results are in complete accord with those obtained by bosonization and by the solution of 
integrable models. 
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I. INTRODUCTION 

A central problcmO in the study of correlated electron 
systems is the classification of the quantum phases and 
critical points of two-dimensional quantum antiferromag- 
nets at zero temperature (T). Of particular interest, be- 
cause of its association with the physics of the cuprate 
superconductors, is the spin S = 1/2 antiferromagnet on 
the square lattice: 



H = J^Si 

(ij) 



(1.1) 



Here Sj are S = 1/2 quantum spin operators on the sites, 
j, of a square lattice, the sum (ij) is over nearest neigh- 
bor links, and J > is the antiferromagnctic exchange 
energy. The ellipsis represents smaller further neighbor 
or ring exchange terms that can be added to deform the 
purely nearest neighbor model: we allow a wide class of 
such terms. Without these additional terms, the ground 
state of H is well known: it is the Neel state with a 
spontaneous staggered magnetic moment 



r]jN e 



T = 0, 



(1.2) 



where rjj = ±1 identifies the square sublattice of site j, 
e is a unit vector with an arbitrary, fixed, orientation 
in spin space, and Nq > is the magnitude of the mo- 
ment. We will refer to this state as the Heisenberg-Neel 
(HN) state, the qualifier denoting the SU(2) (Heisenberg) 
symmetry of H (this allows us to distinguish from cases 
in which the SU(2) spin symmetry has been reduced to 
U(l), which we consider later in the paper). The low en- 
ergy, elementary excitations of the HN state are a doublet 
of gapless spin wave modes, associated with slow modu- 
lations in the orientation e. 

This paper will address the following question and re- 
lated issuesa: what are the possible ground states of H 
which are "near" the HN state ? More precisely, imagine 
slowly reducing the value of Nq (while (1.2) maintains its 
spatial structure) by modifying some of the unspecified 
non-nearest-neighbor terms in H; at some critical point 
in parameter space, Nq will vanish and spin rotation in- 
variance will be restored; how does the ground state of 
H evolve towards and beyond this critical point ? What 
is the field theoretic description of the quantum critical 
point at which Nq vanishes, and of any additional quan- 
tum critical points that may be encountered ? How do 
these results get modified when exchange anisotropies in 
H (like an additional JC^2{ij) Si Z Sj Z term) reduce the 
global spin symmetry of H from SU(2) to U(l) ? 

Along some lines in parameter space, the spatial form 
of (1.2) may change before spin rotation invariance is re- 
stored e.g. the spontaneous magnetization can become 
incommensurate and non-collinear. This is a situation 
of fundamental importance, but in the interests of sim- 
plicity we will postpone discussion of phases associated 
with non-collinear spin correlations to Section VI, where 



(a) 



(b) 



FIG. 1. Sketch of the two simplest possible states 
with bond-centered charge (BC) order: (a) the columnar 
spin-Peierls states, and (b) plaquette state. The different val- 
ues of the Qij on the links are encoded by the different line 
styles; the state in (a) has 3 distinct values of Qy, while that 
in (b) has 2. Both states are 4-fold degenerate; an 8-fold de- 
generate state, with superposition of the above orders, also 
appears as a possible ground state of the generalized action 
for BC order in Appendix |^ 



we present in ( |6.4| ) our 'global' theory for S — 1/2 quan- 
tum antifcrromagnets in two dimensions. The phases and 
phase transitions in Section VI are generalizable to other 



lattices with non-collinear spin correlations e.g. the tri- 
angular lattice. Also, although we will only explicitly 
consider the square lattice antiferromagnet in this pa- 
per, all our results can also be applied to other two- 
dimensional, bipartite lattices, like the honeycomb lat- 
tice. 

Apart from the HN state, states which will play a 
central role in our considerations are those with bond- 
centered .charge order (BC order); these states have been 
arguedBia to be generically contiguous to the HN state, 
and we will also find this here. See Refs. |jM and refer- 
ences therein, for recent numerical evidence for BC order 



in systems like (1.1). The BC states are characterized 
by spontaneous spatial modulations in the expectation 
values of the bond fields 



Qi 



S 7 ; • So 



(1.3) 



where i,j are on nearest neighbor sites. Modulations 
in Qij will lead to corresponding modulations in other 
spin-singlet, charge-conserving observables, including the 
charge density associated with the orbitals that reside on 
the link (in the cuprates, this is the location of the O 
ions). In the simplest cases, the modulation doubles the 
unit cell, with the largest values of (Qij) residing either in 
parallel columns (also referred to as a spin-Peierls state) 
or in disjoint plaquettes (see Fig |l]). 

This paper will address the nature of the quantum crit- 
ical points the system must encounter in evolving from 
the HN to a BC state. On the basis of numerical anal- 
ysis of the square lattice antiferromagnet with first and 
second neighbor exchange interactions, Sushkov, Oitmaa, 
and ZhcngD have proposed that this evolution occurs via 
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an intermediate phase with co-existence of the two order 
parameters i.e. there is first a dimerization transition 
from the HN state to a state with co-existing Neel and 
bond-charge orders (the HN+BC state; see Table | for 
our convention for identifying ground states), and then a 
second transition to the BC state where spin rotation in- 
variance is restored. (Such a scenario was also mentioned 
in the early work of Ref. ||.) As is discussed at length be- 
low, we shall find a number of situations in which our 
field theories display precisely this sequence of transi- 
tions. The other possibility, which general field theoretic 
considerations can never rule out, is a direct first-order 
transition from the HN to the BC state. We will not find 
any clear-cut scenario for a generic (i.e. not multicritical) 
second order transition between the HN and BC states, 
but cannot definitively rule out the possibility that such 
a transition exists (furth er discussion of this point is at 
the end of Section |VC§ ) . 

Another class of our main results concerns the issue 
of spin fra ctio nalization. These will be reviewed below 
in Section [ A . We will find that paramagnetic fraction- 
alized states (F in the notation of Table |) are not con- 
nected by a second-order phase boundary to the HN state 
in the SU(2) symmetric H. In contrast, in spin systems 
with a reduced, easy-plane U(l) symmetry, there carubp-a 
second-order transition between the XY and F statesot!! 

We introduce the computations of this paper by re- 
calling some essential facts about quantum and thermal 
fluctuations near the HN state. At long wavelengths, it is 
accepted that these are well described-by the 0(3) non- 
linear sigma model in 2+1 dimensions!!!!, with the action 

Sa = ^ll [ (5rIl)2 + c2(V * n)2 ] ' (L4) 

where x is the continuum spatial co-ordinate, r is imagi- 
nary time, c is the spin- wave velocity, and g is a (dimen- 
sioned) coupling constant which controls the strength of 
quantum fluctuations. The field n(x,T) satisfies n 2 = 1 
everywhere in spacetime, and denotes the local orienta- 
tion, e, of the short-range Neel order; so 



(1.5) 



However, as the value of N decreases with increasing 7j, 
the action S a ceases to be a complete description of the 
underlying antiferromagnet H . The essential missing in- 
gredients are residual Berry phases associated with the 
precession of the lattice spinsLL3. These appear in the co- 
herent state path integral representation of the quantum 
spin Hamiltonian H in ( p. . l| ) ; they do not have an ob- 
vious continuum limit, and to specify them, we have to 
return to a lattice regularization of S a and its degrees of 
freedom. 

For the spatial co-ordinates, we return to the original 
square lattice associated with H, while we also discretize 
the imaginary time direction: consequently, spacetime 
is replaced by a cubic lattice, on which reside the r\j 
degrees of freedom (the label j is now a three-dimensional 



Label 



HN 



XY 



BC 



Property 



The Hamiltonian is invariant under SU(2) spin ro- 
tations, and spin rotation symmetry is broken by 
(1.2) with No > and e an arbitrary unit vec- 
tor in spin space. There are two gapless spin-wave 
excitations. 



The Hamiltonian has only a U(l) spin rotation 
symmetry about the z axis (d ue t o the presence 
of exchange anisotropics), and (1.2) is obeyed with 



No > and e a unit vector in the x-y plane. There 
is one gapless spin-wave excitation. 



The Hamiltonian has only a U(l) spin rotation 
symmetry about the z axis (d ue t o the presence 
of exchange anisotropics), and (1.2) is obeyed with 
No > an Ising order parameter, and e a unit vec- 
tor along the ±z directions. There are no gapless 
spin-wave excitations. 



Lattice space group symmetry is broken by a spon- 
taneous modulation in the values of the bond- 
centered charge order parameter defined in 
(1.3). The simplest patterns, in Fig 0, double the 
unit cell to two sites. 



The ground state has topological order (leading to 
a four-fold degeneracy on a large torus) and frac- 
tionalization. There are one or more gapped exci- 
tations with spin S z = ±1/2, and also a spinless 
Z2 vortex excitation. 



Spin rotation symmetry is broken by a non- 
collinear, polarization of the spins. There are three 
gapless spin-wave excitations if the Hamiltonian 
has full SU(2) symmetry, and only one if the sym- 
metry is U(l). 



TABLE I. We identify ground states of two-dimensional 
antiferromagnets by a label which specifies all of the proper- 
ties from the above list which are satisfied by the state. It is 
assumed that if a property not contained in a state's label, 
that state does not satisfy its requirements; so e.g. while both 
the HN and HN+F states obey (R~ 



1.2), the HN state does not 



obey property F, while the HN+F state does. Along the same 
lines, the F state is invariant under spin rotations, as it does 
not obey properties HN, XY, I, or P. 
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cub ic la ttice index, rather than a two-dimensional index 
in (LI) — the interpretation of j should usually be clear 
from the context). On this cubic lattice, the Berry phase 
term in the partition function is e~ SB with 



i r/jAjr, 

j 



(1.6) 



where, as before, rjj = 1 on one of the square sublattices 
and rjj = — 1 on the other (note rjj is independent of r), 
and the leading i = The quantity Aj^ is defined 

to be half (for S = 1/2) the signed area of the spherical 
triangle formed between ny, and no, with no an 

arbitrary reference unit vector (the index fi extends over 
the spacetime directions x,y,r). We can choose n at 
our convenience, and a convenient choice is usually no = 
(0, 0, 1). The explicit expression for half the area, A, of 
the spherical triangle bounded by 111,2,3 ial3 

lA _ 1 + ni • n 2 + n 2 • n 3 + n 3 • ni + mi • (n 2 x n 3 ) 
{2(1 + m • n 2 )(l + n 2 • n 3 )(l + n 3 • m)} 1/2 

(1.7) 

It can be shown from this definition, or more geomet- 
rically from the interpretation of A as a spherical area, 
that changes in the choice of no amount to a gauge trans- 
formation of the Ajp, with 



A 



.hi 



A 



hi 



•jj+ii 



"3> 



(1.8) 



where <j)j is half the area of the spherical triangle formed 
by rij and the old and new choi ces f or no . It is easily seen 



that (L6) is invariant under (L8) provided we choose 
periodic boundary conditions in the r direction — this we 
will always do. As a companion to Sb, we also write 
down the continuum S a in ( |l.4| ) on the same hypercubic 
lattice: 



(1.9) 



where the fi are vectors of length equal to the lattice spac- 
ing, and extending over the x, y, and r directions. The 
dimensionless coupling constant g is proportional to g 1 
while c has been absorbed by the choice of the temporal 
lattice spacing. For the above discretizations to be mean- 
ingful representations of the underlying antiferromagnet, 
we exclude values of g so large that the n correlation 
length is of order a lattice spacing or smaller. 

In the analogous formulation of one-dimensional 
antiferromagnetal3, a natural continuum representation 
of the Berry phase term does exist. A careful examina- 
tion of Sb for smooth configurations of n,, shows that it 
can be written as 



id f , , / dn dn 
— / dxdm ■ — x — 

47T / \ OT OX 



(1.10) 



with 9 = 7r. This is a topological ter m, an d the contin- 
uum theory S a + Sg (the integral in (1.4) now extends 



over one spatial dimension) was used by Haldane to ar- 
gue for the fundamental distinction between integer and 
half-integer spin chains. While the continuum theory has 
been valuable in establishing this distinction, very little 
information on the physics at 9 = n has been obtained di- 
rectly from the continuum functional integral of S a + Sq. 
Rather, the most efficient approach has been to refer back 
to the lattice S=l/2 spin chains, and to analyze them by 
the methods of abelian bosonization. This paper will 
perform a direct analysis of the lattice actions 5* n + Sb, 
in both spatial dimensions d — 1 and d = 2. We will use 
parallel methods in the two cases, and the comparison 
with the known bosonization results in d = 1 will act as 
a significant consistency check of our conclusions. 

Although it is not essential, it is convenient to intro- 
duce another term in the final form lattice action we shall 
study — this has the advantage of yielding an independent 
dimensionless coupling constant (in addition to g) which 
can be used to explore a potentially richer phase diagram. 
We wish to write down a term associated with the gauge 
potential Aj^. Any such term m ust c learly be invariant 
under the gauge transformation (|l.8|) . Further, as only 
the values of n., are observable, ( |1.7[ ) indicates that the 
term should be periodic under Aj^ — > Aj^ + 2ir. These 
requirements strongly constrain the allowed terms, and 
the simplest permissible one is 



S'a 



P 2 



(1.11) 



This term is written in notation standard in lattice gauge 
theory: the sum over □ extends over all plaquettes of the 
cubic lattice, the indices jj,, z/, A extend over the x, y, r 
directions. The symbol A M represents a discrete lattice 
derivative (A M /j = fj+p, — fj), and so Sa depends on 
the A flux threading each plaquette. In principle, we 
should allow for different coupling constants associated 
with fluxes in the spatial and temporal directions, but 
we have chosen them equal for convenience. As we will 
discuss shortly, an importan t too l in our analysis is a 
dual lattice representation of ( 1.11 ) . It is somewhat more 
co nvenie nt to perform this in a "Villain representation" 
of ( 1.11 ), which replaces the exponential of a cosine by a 
periodic sum over gaussians. So we shall replace S' A by 

Sa = ^2 I>^A A ^A - 27r ^) 2 - ( L12 ) 

where is an integer-valued vector field on the sites, 
J, of the dual cubic lattice; q couples to the A flux on 
the plaquette it pierces. It is also clear that SU becomes 
singular as e 2 — > 0, and we shall therefore exclude this 
limit from our considerations; we are primarily interested 
in moderately large values of e 2 , where the connection 
to Heisenberg antiferromagnets is evident. However, the 
limit e 2 — > oo is not expected to be singular, and the 
properties at e 2 = oo should be smoothly connected to 
those for large e 2 . 
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HN 



HN+BC 



BC 



g 



1c 



g 2c 



g 



FIG. 2. Proposed phase diagram of the model Z in ( LIS ) 
as a function of g for a range of e 2 . The transition at g = gi c 
involves onset of BC order like that in Fig |l[ this is in the uni- 
versality class of the 2+1 dimensional Z4 clock model, which 
is in turn described by the <p 4 field theory with 0(2) symme- 
try. The HN order vanishes at g — g2 C > gic with a transition 
described by the 2+1 dimensional ip 4 field theory with 0(3) 
symmetry. A direct first order transition between the HN and 
BC phases is also possible. We have not find any clear-cut 
scenario for a direct second order transition between the HN 
and BC states, but cannot definitively rule out the possibility 
that such a transition exists for some values of e 2 (further 



discussion of this point is at the end of Section V C \\ ) 



It is useful to collect all the terms introduced so far. 
One of the central purposes of this paper is to understand 
the phase diagram of the partition function 

z = J II dn Atf - !) CX P (Sn -s A - S B ) 



{«»} j 



(1.13) 



in the plane of the coupling constants g and e 2 . Z consti- 
tutes a "minimal model" for a theory of S = 1/2 quantum 
antiferromagnets in two dimensions: it incorporates the 
essential ingredients — quantum fluctuations of the short- 
range Neel order represented by the non-linear sigma 
model plus the Berry phases accounting for the quan- 
tization of spin. This model, and the duality methods 
we shall employ, are closely related to those introduced 
in Ref. [I| see also Appendix [A| where we introduce an 
alternative formulation of Z using fields in CP 1 — this 
formulation is more convenient for some purposes, and 
is closer to the approach of Ref. ^. The presence of the 
complex Berry phases associated with Sb means that a 
Monte Carlo evaluation of Z will not be straightforward; 
nevertheless, we hope that such evaluations will be at- 
tempted by experts in the future. Our strategy will be to 
employ a variety of duality and field-theoretic techniques 
to delineate the phase diagram of Z. We will present 
strong arguments that only phases with HN and BC or- 
der appear in this phase diagram, and that there are no 
states with spin frac tion alization in the SU(2) symmetric 
Z. (see also Section I A below). A likely scenario for the 
sequence of phases as a function of increasing g is HN, 
HN+BC, BC, for all non-zero values of e 2 , as shown in 

As we have already indicated, the connection of our 
models with the successful theory of one-dimensional an- 
tiferromagnets will be an important guide to our analy- 
sis. For these systems, a complete understanding of the 
physics was achieved by considering the phase diagram 
in a wider space: the spin-symmetry was reduced to U (1) 
by introducing an easy-plane or easy-axis anisotropy, 



and the phase diagram was studied as a function of the 
strength of the anisotropy. After mapping to continuum 
field theories by abelian bosonization methods applica- 
ble for systems with only a U(l) symmetry, it was found 
that one could nevertheless identify lines in the parame- 
ter space of the field theories where the Heisenberg SU(2) 
symmetry was restored. Thus one achieved an under- 
standing of SU(2) symmetric antiferromagnets by com- 
pletely exploring the phase diagram of less-symmetric, 
and simpler, systems. A similar strategy will be an im- 
portant ingredient in our arsenal here: we will reduce 
the symmetry in the n space to U(l), by introducing ex- 
change anisotropy in S n and so d eforming Z to ^u(i) 
(the explicit form is below in (1.15)). We find then that 
duality transforms of -Ztj(i) can be carried through to 
completion. This yields models which we will study by 
mappings to continuum field theories in 2+1 dimensions 
which are amenable to standard renormalization group 
analyses. We will then search for lines in the parame- 
ter space of these dual models at which SU(2) symmetry 
is present; the properties of such lines will allow us to 
reach important conclusions about the SU(2) symmetric 
partition function Z. 



A. Fractionalization 

Much attention has focused on ground states of two- 
dimensional antiferromagnets with an "order" quite dis- 
tinct from those discussed so far. These are states with 
fractionalized spin excitations (spin S = 1/2) and associ- 
ated topological, spin-singlet, excitations. States with 
this pfGipeirty are denoted by F, see Table |). In a 
theoryliaila of paramagnetic F states in systems with 
SU(2) symmetry, fractionalization was linked with non- 
collinear spin correlationsliH In particular, restoring 
spin rotation invariance in a state with long-range, non- 
collinear, magnetic order (a P state, see Table m by a 
second-order quantum phase transition led naturally to 
a paramagnetic F state. Suc h P s tates are not expected 
to appear in the model Z in ( |l.l3| ), and so the theory of 
Ref. 14 leads to the conclusion that the SU(2) symmetric 
Z does not contain the F state either. Our results here 
will support this conclusion. 

As already noted above, we will also study a deforma- 
tion of Z to anisotropic spin models, ^u(i) with only a 
U(l) spin symmetry. For easy-plane anisotropy, this will 
deform the HN state into the XY state (see Table |jj) . In 
such a situation, the rij field will fluctuate primarily in 
the x-y plane in spin space, and it is therefore useful to 
parameterize in terms of an angular variable 6j by 



(cos 9j , sin 9j , 0) . 



(1.14) 



We will show in Sections [fVB andjv that in such an 
easy-plane limit, the model Z in (1.13) is equivalent to 
the following XY model coupled to a Z2 gauge theory in 
2+1 dimensions: 
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z u(1) = ]r f n ^ cx P ( ^ 5^ n sjj+h 



.9 



E" ( L15 ) 

i / 



with the coupling K is related to e 2 by (4.11); actually 
the precise equivalence is to a model in which the cosine 
in (|1.15|) is replaced a periodic Gaussian- Villain form, as 
in ( [4. 10D . The Z<x gauge field sjj+p, resides on the links 
of the direct lattice; apart from the usual Maxwell term, 
it also carries the las t Berry phase term which i s the 
remnant of Sp in (F6). Models closely related to ( 1.15 ) 
have recently, been studied by Senthil and Fisher and 
collaborators! - !!!]: we have thus established a surprising 
connection between their work and the easy-plane limit 
of model Z of Ref. || and the present paper. 

We make a brief aside to place the work of Senthil and 
Fisher in the context of the models being studied here: 
their work has shed new light on subject of fractional- 
ization injspin systems with U(l) symmetry. In their 
first papera, Senthil and Fisher considered a model with 
charge fluctuations and SU(2) spin symmetry, but ex- 
amined fractionalization of the Cooper pairs associated 
with superfluid order and the U(l) charge symmetry. Al- 
ternatively, we could replace the Cooper pair bosons by 
the hard-core bosons associated with S=l/2 spin systems 
(with the spin raising operators S + = ~ e z6> , the oper- 
ator creating the hard-core boson (see ( 4.21[ ) later)), and 
transcribe their results to quantum spin systems with 
only a U(l) symmetry. Precisely such a point of view 
was taken in Section II of Ref. |l(J. 

Our analysis of ^u(i) w nl follow and generalize the 
framework laid by these studies of Z^ gauge theories. 
Duality transformations of the model -^u(i) show that, 
unlike Z, -Zu(i) does indeed contain a paramagnetic F 
state in two dimensions; further, it appears possible to 
directly connect the XY and F states by a second-order 
quantum phase transition. The appearance of these F 
phases is intimately linked to the half-angle variab le ap - 
pearing as the argument of the cosine term in (1.15); 
indeed e %e l 2 is the 'square root' of the S z = 1 spin-flip 
boson &t, and so the former is the creation operator a 
S z = 1/2 spinon. For large K and large g, the fluctua- 
tions of the Zi gauge field 8j j+p and the XY order will 
both be suppressed, and this will introduce a paramag- 
netic F state in whi ch su ch spinons can propagate freely 
by the 1/g term in ( 1.15 ). 

Our main interest here is in spin syst ems with SU(2) 
symmetry, and so in Sections IVB 2 and V C we will con- 
sider a wide class of generalizations of -^u(i) m the hope 
of identifying special regimes of parameter space where 
the full SU(2) symmetry may be restored. Our search 
for such spaces for enhanced symmetry will be success- 
ful: one of our main findings will be that such F states 
do not survive the restoration of SU(2) symmetry. Along 
with other reasons to be discussed in Sections [n] and III , 



thi s lead s to our conclusion that the SU(2) symmetric Z 
in (1.13) does not contain any fractionalized states. Sim- 
ilar reasoning will also allow us to conclude that the HN 
and F states cannot be separated by a single, second- 
order quantum critical point, in contrast to that found 
in ^u(i) between the XY and F states. 

To further clarify the issues of fractionalization in sys- 
tems with SU(2) symmetry, we will introduce a distinct 
lattice model, Zp, in Section VI. This model generalizes 
Z to explicitly allow for a non-collincar, magnetically or- 
dered state (a P state) . We will argue that it also contains 
an F state, and discuss the structure of the phase diagram 
of Zp. This will allow us to delineate the routes by which 
it is possible to connect states with HN, BC, and/or F 
order in systems with full SU(2) spin symmetry. We ex- 
pect the conclusions in this section to be rather general, 
and apply to a wide class of SU(2) symmetric models 
with phases characterized as in Table [j]. 



B. Outline 

This is a lengthy paper, and so it is helpful for the 
reader to have an overview of the logical relationships 
between the sections at the outset. We will begin in 
Section |n| by a direct assau lt on the SU(2) invariant par- 
tition function Z in ( 1.13| ). We will argue that BC or- 
der appears at large g and then proceed in Section III 
to a discussion of how the system interpolates between 
the small g HN phase and the large g BC phase. We 
will present a number of arguments based upon the sta- 
tistical mechanics of defects in the <p 4 field theory in 3 
spacetime dimensions that support the proposal of an 
intermediate phase with co-existing HN and BC order. 
The next two sections take the reader on a long detour. 
In Section IV we turn our attention to one dimensional 



antifcrromagnets and allow for the system to acquire a 
uniaxial anisotropy and a reduced U(l) symmetry. For 
such systems, we are able to develop rather precise du- 
ality methods which completely account for the defects 
and their Berry phases; we are also able to restore SU(2) 
symmetry along special lines in parameter space. The 
results obtained in this manner are found to be in excel- 
lent accord with those obtained earlier by bosonization 
methods. Emboldened by this success, we proceed with 
the application of the same methods in two dimensions 
m Section [v|. The main results again support the co- 
existence of HN and BC order in the SU(2) invariant 
systems, in agreement with Section III. Finally, in Sec- 
tion VI we return to the main road, and focus direct 



attention on SU(2) invariant antifcrromagnets in two di- 
mensions. We introduce a generalization of Z, denoted 
Zp, which explicitly accounts for non-collinear spin cor- 
relations expected in frustrated antiferromagnets. This 
model does contain paramagnetic fractionalized phases, 
and we present phase diagrams which describe competi- 
tion between HN, BC, and F orders. 



G 



II. DUALITY AND BOND ORDER 



The model Z in (1.13) has a clear association with the 
HN state: in the limit of small g, the coupling in S n 
prefers that all the rij align along a common direction 
in spin space, and this leads to a ground state with HN 
order, as in (1.2). This section will introduce a dual- 
ity transformation of Z which demonstrates its intimate 
connection with BC states. We will argue that BC order 
appears for large g in the SU(2)-symmetric Z. We will 
also indicate why there is no such requirement for Zuny 

This first s tep in the duality transformation, as usual, 
is to rewrite (1.12) by the Poisson summation formula: 



E ex p f - 



(2.1) 



where (like q^) is an integer- valued vector field on 
the links of the dual lattice. Here, and henceforth, wc 
drop overall normalization co nsta nts in front of partition 
functions. The last term in ( [2.l| ) has the structure of a 
'Chern-Simons' coupling, and connects aj M to the A flux 
on the plaquette of the direct lattice that it pierces. 
Next, we write Sb in a form more amenable to duality 



transformations, 
which satisfies 



Choose a 'background' aj M = a)j flux 



e Ml/A A„a" A = rjjSfir, (2.2) 

where j is the direct lattice site in the center of the pla- 
quette defined by the curl on the left-hand-side. Any 
integer- valued solution of (2.2) is an acceptable choice 
f or a " , and a convenient choice is shown in Fig ||[ With 
( |2.2| ), the Berry phase term also takes the form of Chern- 
Simons coupling: 



: E 



□ 



J'A 



(2.3) 



We can now combine ( |2.1| ) and ( |2.3| ) and, after a shift of 
a jv by ct° M , obtain a new exact representation of Z: 

Z = J II °MK 2 - !) exp(-Sn - S a 

i ^ tuvxa-on k-v-Ajx j , (2.4) 



where 



J4> 



The structure of this representation of Z also allows us to 
give a simple physical interpretation of aj„ . The quantity 



x 
x 

X 



x 

X 
X 



FIG. 3. Specification of the non-zero values of the fixed 
field dj M . The circles are the sites of the direct lattice, j, 
while the crosses are the sites of the dual lattice, j; the latter 
are also offset by half a lattice spacing in the direction out 
of the paper (the fj, = t direction) . The all zero for 

fi = t,x, while the only non-zero v alue s of shown 
above. Notice that the a flux obeys (p.2[). 



( a JM — a j/i) 2 f° r I 1 = x >y ^ s a term in action associated 
with a spatial link on the dual lattice, and is therefore 
surely associated with a measure of the energy of the 
bond on the direct lattice which intersects it. We may 



therefore identify the bond variable in (1.3) as 



(a Jy - a° Jy ) 2 , 



(2.6) 



where the direct lattice link on the left-hand-side inter- 
sects the dual lattice link on the right-hand-side. A sim- 
ilar relation holds for x «-> y. 

We have now assembled the ingredients necessary to 
address the physics of the large g regime. Here the n 
field will fluctuate strongly, and so it appears natural to 
integrate it out. To do this, we need to evaluate the 
average 



W{a, 



exp ^—i e^xajn^-uAjxj 



In terms of lF(a M ), the partition function is 



z = w m ex p 



(2.7) 



(2.8) 



Strong fluctuations in n will lead to concomitant fluctu- 
ations in the flux e^xA^Ajx over all real val ues between 
and 2-7T. This means that the average in (2/7) will be 
strongly dominated by configurations in aj M which obey 



^tiuX^vO-jX — 0. 



(2.9) 



(2.5) In other words, the large g region is in a "a-Meissner" 
phase which expels the flux of ajM. We will restrict at- 



tention to Oj M configurations which satisfy (2.9) in most 
of the remainder of t his section. A more careful discus- 
sion of the origin of (2.9) is given in Appendix H; there 
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we show that the configurations which violate ( |2.9| ) lead 
mainly to a renormalization of the coupling constant e in 
S a . The discussion in Appendix |b] also shows that the 
analysis of this section also applies to models which a 
bare value e = oo; the renormalized e appearing in S a is 
finite even in that case. 

This is also a convenient point to make a brief aside 
on systems with an easy-plane anisotropy, and a reduced 
U(l ) sym metry. Here n lies in the x-y plane as indicated 
in ( 1.14 ), and the flux e^xA^Ajx on any plaquette is 
typically zero. Only if plaquette happens to contain a 
vortex in the angle 9 is the flux non-zero, and then it 
takes the discrete values ±7r (half the area of a hemi- 
sphere). More generally, allowing for multiple vortices, 
the allowed values of the flux e^xA^Ajx are mir with m 
integer. In the large g region, we can assume that there 
were be strong fluctuations in the value of m in the pla- 
quettes, and so the average in ( |2.7| ) will be dominated by 
configurations in which satisfy 



e^xAvdjx = (mod 2), 



(2.10) 



and the flux of a ^, is only expe lled modulo 2. The dif- 



ference between (2J)) and ( 2.10 ) is the key reason from 
the strong distinction between the properties of SU(2) 
and U(l) symmetric spin systems. We will show in Sec- 



tion VA that the partitio n func tion (2.8), when evalu- 
ated under the constraint ( 2.10| ), reduces exactly to an 
Ising model on the cubic lattice in D = 3 dimensions, 
with the signs of the couplings chosen so that every spa- 
tial plaquette is frustrated; this Ising model is dual to 
the Ising gauge theory plus Berry phases obtained in the 
large g limit of ( 1.15| ). Precisely the same Ising model 
was considered in Refs. ^6|JT^ in a context somewhat dif- 
ferent from the present, but closely related to the models 
of frustrated antiferromagnets to be considered in Sec- 
tion 0. We can also already notice here that double 
vortices (with m = ±2) don't couple at all to fluctua- 
tions in aj M , and so can selectively proliferate for large g: 
this is the effect responsible for appearance of F states in 

Returning to SU(2)-invariant systems, and the a- 
Meissner phase obeying (2J)), we can solve the constraint 
by writing 



(2.11) 



where hj is a "height" on the sites of the dual lattice. 
The partition function Z is now equivalent to the three- 
dimensional height-model partition function 




e 2 



3 



(2.12) 



The model 
of occasion 



been studied previously on a number 
and we now state the main results. 
The most important property of Zh is that the heights 
are in a smooth phase for all values of e; in other words, 



any state of Zh has a definite value for its average height 
(hj) , where the average is over both quantum fluctuations 
and over sites of the dual lattice — this is a generic prop- 
erty of height models in three dimensions. Further more , 
the combination of the background a?^ terms in (2.12) 
and the existence of an average height implies that the 
translational s ymm etry will be broken by a modulation 
in the Qij in ( |2.6|) . In other words, the ground state of 
Zh has BC order. 

A review of the prior analyses of Zh which led to the 
above results is provided in Appendix There, we also 
establish the close connection between Zh and a num- 
ber of other models for quantum paramagnets that have 
appeared in the literature. In particular, there is a sim- 
ple and direct|-Gppaection between Zh and the "quantum 
dimer" modeled E3. Furthermore, Zh can also be writ- 
ten as a Coulomb gas of point instanton charges, each of 
which carries a Berry phase: these instantons are noth- 
ing but 'hedgehog' defects in the n field, and the Berry 
phases implied by Zh are precisely equal to the hedgehog 
Berry phases computed by HaldanoE£l. Finally, Zh can 
also be written as a frustrated sine-Gordon model. This 
last form is the most useful in determining the ground 
states of Zh- This simplest state which emerges is the 
4-fold degenerate columnar BC state shown in Fig [j]a. 
However, within the generalized— earameter space of the 
sine-Gordon model we also findc^l the 4-fold degenerate 
plaquette ground state of Fig |l|b, and an 8-fold degener- 
ate states involving co-existence of the orders in Fig |l|a 
and b. The symbol BC refers collectively to any one of 
these states; an antiferromagnet could also have quan- 
tum phase transitions between different BC states-these 
transitions can also be discussed easily within the frame- 
work of the frustrated sine-Gordon modelcJ, but we will 
not dwell on them here. 



III. CONNECTING THE NEEL AND BOND 
ORDER STATES 

We have so far exami ned t he phases of the SU(2) sym- 
metric model Z (Eqn ( 1.15 )) in the limits of small and 
large g. For small g we have the HN state with broken 
spin rotation symmetry, while for large g we have SU(2) 
symmetric ground states with BC order. This section will 
comment on the evolution of the ground state as g moves 
between these two limits. Our anal ysis will focus of the 
behavior of the function W(a^) in (2.7) as a function of 



g: this will permit us to study the influence of the HN 
order parameter n on the fluctuations of the q^„, w hich 



control the BC order. We will begin in Section HI A with 
limit of small g, where n fluctuations can be co mpute d in 



the spin- wave approximation. Next, in Section III B 



we 



will consider the point where n fluctuations are critical 
i.e. at the point where SU(2) symmetry is broken, and 
describe the associated structure of W(a tl ). The criti- 
cal point is a strongly coupled theory in three spacetime 
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dimensions, and so our arguments here do have ad hoc 
steps which are based mainly on physical arguments. 



A. Spin wave expansion 

First, we describe the small g spin-wave expansion. 
Here, we assume (n) = A*o(0,0, 1). We perform the fa- 
miliar spin- wave expansion of the non- linear sigma model 

by parameterizing n = (tt x , n y , ^Jl — n x — n y ) and ex- 
panding the action, including the expressions for Aj^ in 
powers of ir x ,ir y . Integrating out the tt x , Tr y fluctuations 
to lowest order in g we obtain 

W(a^,) = exp(-£ stl) (a M )) 
9 



[G(i- j )G(i-j + fi-P)-G(i-j + pL)G(i-j- 
d 3 k ge lk -i/2 



G(j) 



8ir 3 3 — cos k x — cos k y — cos k z 



(3.1) 



We now need to study the fluctuations, as controlled 
by the action S a + S sw , to determine the fate of the bond 
order in this framework. A closely related action was 
numerically simulated earlier in Ref. with an action 
in which the double summation in S sw over i and j was 
replaced with a single summation over the on-site term 
with i = j, and p = v. Numerical evaluation of the 



expression for G in (3.1) shows that this on-site term is 
over 10 times larger than the off-site terms, and so the 
truncation involved in mapping to the model of Ref. |^ is 
reasonable. The primary res ult o f Ref. ||[ which we expect 
to apply to the full action in (3.1), was that there is a crit- 
ical value of g above which there is an onset of BC order. 
However, the present small g expansion does not allow us 
to determine if this onset occurs in a regime where the 
small g spin-wave expansion about the HN state is valid. 
If the transition in S a + S sw occurs in its regime of appli- 
cability, then the ground state has HN+BC order beyond 
the critical point. So the structure of the spin- wave ex- 
pansion in powers of g is compatible with both the HN 
and HN+BC states, but does not allow us to definitively 
conclude that the HN+BC state must exist. 



B. Critical ip 4 field theory 

To study the competition between the HN and BC 
states further, we imagine turning up the value of g to 
the critical point where SU(2) symmetry is first restored, 
and No vanishes. This section will make the assumption 
that the fluctuations of n at this critical point are de- 
scribed by the field theory associated with the critical 
point in S n itself i.e. the 3-component ip 4 field theory in 
2+1 dimensions with 0(3) symmetry. We will argue, by 



theory, that this assumption is justified only if HN order 
vanishes at a transition between the HN+BC and BC 
states. In other words, we will find that a critical point 
described by the 3-component p> 4 field theory has finite, 
non-critical BC order. The alert reader will notice that 
these arguments do not logically exclude the possibility 
that there is some other unknown critical field theory de- 
scribing the vanishing of HN order, and which describes 
a second-order phase transition between the HN and BC 
states; however, we will not find any likely candidates for 
such a critical theory here. 

To evaluate W(a fJi ) we need some understanding of the 
fluctuations of the A flux. As we saw in the large g theory 
in Section [0], these are controlled by the fluctuations of 
the hedgehog point-defects in the n fluctuations. So as 
a first step, we should determine the correlations of the 
hedgehogs at the critical point of the ip 4 field theory. 
We will be able to make a number of statements about 
these correlations in an expansion in e = 4 — D (where 
D is the dimensionality of spacetime) by an extension of 
methods developed by HalperinEj some time ago, which 
are described in Appendix 19. We begin by outlining the 
physical ideas behind our analysis, and will then apply 
them to the computation of W(a M ). 

It is useful to first contrast the critical behavior of the 
hedgehogs with a well-known and familiar example: vor- 
tices in the D = 2 XY model which exhibit the Kosterlitz- 
Thouless (KT) transition. Imagine we are examining a 
model with a short-distance momentum cutoff, A. Let 
the mean density of vortices be p v , and a dimensionless 
measure of this is p v A~ 2 . Now we perform an RG trans- 
formation, integrating out tightly-bound dipoles of vor- 
tices, and gradually reduce the value of A. It is known 
that at the KT critical point the vortex fugacity ulti- 
mately flows to zero; in other words, as we scale A to 
smaller values, the dimensionless vortex density, p v A~ 2 , 
ultimately scales to zero. This is the precise form of the 
statement that there are no free vortices at the KT crit- 
ical point. The behavior of the density of hedgehogs, 
Ph, at—the. critical point of the D = 3 0(3) model is 
knowr£j'Ej to be dramatically different: p~hk~ 3 remains 
a finite number of order unity even after RG scaling to 
the longest scales. So, loosely speaking, there is finite 
density of "free" hedgehogs at the critical point. We 
will compute correlations of this critical ensemble of de- 
fects of below. Our main physical observation will be 
that these c ritic al de cay correlations decay sufficiently 
rapidly (see (3J3) and (3/7) below) that there is little ma- 
terial difference from the corresponding correlations in 
the paramagnetic phase. The consequences of the latter 
correlations were explored earlier in the large g theory: 
we saw in Section || that they led to BC order in the 
paramagnetic phase. So it is reasonable to conclude that 
finite BC order is also present at the critical point of the 
HN order. 

Let us now describe the correlations of the hedgehogs 



at the critical point of the 0(3) ip field theory in 2+1 
an evaluation of W(a ll ) at the critical point of the ip 4 field dimensions. Let Ph(r) be the topological charge density 







operator of hedgehogs at the spacetime point r (so ph 
(\ph{r)\))- We consider first the two-point correlator 



C h (r) = ( Ph (r) Ph (Q)) 



(3.2) 



The topological character of the Ph(r), and its association 
with a continuum field ip{r) implies that this correlator 
obeys an overall charge neutrality condition 



d 3 rC h (r) = 0. 



(3.3) 



We explicitly evaluate Ch (f) in an expansion in e = 4 — 
D in Appendix The topological character is special 
to D — 3, and so the constraints of the conservation 
law (3.3) are not reflected in this computation. We find 



that the long-distance decay of Ch (r) is controlled by the 
scaling dimension 



dim[p h (r)] = (9 + %)/2 
where r\h has the e expansion 
r] h = 3r] + 0(e 3 ) 
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0(e 3 ), 



(3.4) 



(3.5) 



where rj is the anomalous dimension of the field ip. This 
scaling dimension implies that Ch has the long distance 
behavior 



C h (r — > oo) 



1 



| r |9+% ' 



(3.6) 



One might question the applicability of such an expan- 
sion to D = 3, given the neglect of the conserv ation law 
(3.3) to all orders in the e expansion. Indeed, if ( |3.3[ ) was 
obeyed by the scaling limit of the field theory, then we 
would expect dim[p/j(r)] = 3, or % = —3. We think this 
is extremely unlikely, and the more plausible scenario is 
that the conservation (3.3) is not obeyed by the scaling 
limit field theory (this is analogous to the the appearance 
of an "anomaly" in a field theory): to properly define 
the field theory, some consistent short-distance regular- 
ization has to b e introduced, and this regularized theory 
does obey (3.3) after including short distance contribu- 
tions which are beyond the scaling limit. The reason 
for the presence of the "anomaly" becomes clear after 
implementing the consequences of (3.3) and (3.6) in mo- 
mentum space: Ch(k), the Fourier transform of C^(r), 
has the following small k expansion 



C h (k -> 0) ~ Cl k 2 + ... + ca* 8+n * + 



(3.7) 



where the first set of ellipsis in (f3J) refer to various an- 
alytic terms in integer powers k 2 , and Ci ( 2 are constants. 
Even for rjh = —3, note that the second singular term 
remains subdominant to the analytic term proportional 
to c%. The latter is controlled by non- universal short- 
distance effects, and it appears that these are the pri- 
mary degrees of freedom controlling the conservation of 



total topological charge. For the scaling limit theory to 
be tightly constrained by the conservation law, we would 
require c\ — 0, something that does not seem likely in 
general, given its dependence on non-universal features. 
So we conclude that the critical theory can ignore the 
conservation law, and the e expansion for r\h remains ap- 
plicable in D = 3. The singular piece of Ch decays very 
rapidly with r, and the small momentum behavior of Ch 
is dominated by a non-universal, analytic term with co- 
efficient ci =/= 0: these are the primary conclusions on 
the properties of Ch(r) that we will need in our analysis 
below. 

The e expansion also allows computation of higher or- 
der correlations of ph- To leading order in e and at small 
k, we can neglect higher order cumulants in the correla- 
tions of ph- So we will assume that the fluctuations of ph 
are controlled by the distribution controlled only by the 
second cumulant, which is 



exp 



Mfc)l 2 

C h (k) 



(3.8) 



It now remains to evaluate the distribution of the A 
flux associated with the above distribution of ph, and to 
thence obtain W{a ll ). Let us define the 'electromagnetic' 
field 



(3.9) 



As we are focusing attention on long-wavelengths here, 
we have replaced the lattice derivative by a continuum 
partial derivative. In the presence of a fixed hedgehog 
density ph, the £ field obeys the 'Maxwell' equations of 
motion 







(3.10) 



From these equations, and from ( |3.8| ), we can deduce that 
the distribution of £ is controlled by 



exp 



It. 



k 2 \£,{k)\ 2 
47r 2 C h (fc) 



(3.11) 



At this point, it is temptin g to compute W(a 



from the distribution (3.11) and to conclude that 



exp 



4e 



4ir 2 C h (k) 



(3.12) 



However, this expression is not invariant under the 
integer-valued gauge transformations aj M — * — AJj 
which are respected by the original definition in (|2 7|) . 
To remedy this, we will have to put ( [3.11 ) back on the 
lattice, and compute the average over £ M while respect- 



ing the periodicity Aj 



2tt. This we do in the 



following paragraph. Readers not interested in this tech- 
nical step may skip ahead to the main result ( p. 15 ) which 
replaces (3.12). 
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First we decouple the distribution in (3.11) by a 



Hubbard-Stratonovich field 
(VS, exp L\ £ l^iM|5,(fc)| 2 + C£ ^{k)S,{- 

J \ k k 

(3.13) 

Then we restore t he pr oper flux-periodicity in the last 
term by replacing (3.13) by 



E I ^S^exp 



^ 2 C h {k) 



2^ 

k 



sites of the dual lattice, 
result, 



E ex p -E 



27T 2 C h (k) 
k 2 



a^k)- A M /i(fc)| s 



for a M becomes precisely that in (B4) in Appendix 



the 0(3) ip 4 model; the BC order is finite on both sides 
of this transition, as well as the critical point, and does 
not play an essential role in its critical theory at g 2c ■ 



k) 



(3.14) 



We can now e valua te W(a tl ) with the distribution of A 



controlled by ( 3.14 ): the sum over qj^ implies that <Sj M 
must be an integer, while the integral over A implies that 
<%jii = cijp. — A M /ij M where h is integer- valued field on the 



From this we obtain our main 



(3.15) 



for W{a^) at the criti cal point of the 0(3) ip 4 field theory. 
The expression (3.15) modifies ( 3.1 2j ) to allow invariance 
under integer- valued gauge transformations of the a^. 

The status of the BC order at this critical point is now 
controlled by the a M fluctuations associated with the dis- 
tributi on e~ 5a W(a fJ/ ). It is surely acceptable to replace 
Ch in fl3.15|) by its leading small k expansion in (3.7): 
Ch(k) = cifc 2 , with ci ^ 0. Then, the distribution 



the case of the large g limit. Our conclusions here are 
therefore the same as those in the discussion following 
(B4): the systems expels a flux, and is in the a-Meissner 
phase. There is, therefore, well-developed BC order at the 
critical point of the 0(3) ip 4 model. The BC order fluctu- 
ations are also not critical and this justifies, a posteriori, 
that the fluctuations of the HN order at the critical point 
at which No vanishes are described by the 0(3) ip 4 model 
alone. 

This section has therefore dilineated an explicit route 
for the destruction of the HN order by second-order quan- 
tum transitions, which was sketched in Fig g. With in- 
creasing g, and before the HN order vanishes, there is an 
onset of BC order at g = g± c ; for BC order of the two 
types shown in Fig [|, this transition is described by the 
i?4 clock model, which is in turn in the universality class 
of the 0(2) ip 4 field theory. For g > gi c , the system is 
in the HN+BC state. At a larger value of g = g 2c HN 
order disappears by a continuous transition described by 



' 1. U(l) symmetry 

We mention in this subsection the application of the 
analysis above to systems with only a U(l) spin sym- 
metry. We will consider such models in much detail in 
Section ^ using rather more precise methods. So we can 
test our present approach and see if consistent conclu- 
sions are obtained. For Z\jn), as we already noted in 
Section y, the 'flux' £^ = e^xA^Ax through any pla- 
quette is simply n times the vortex number piercing that 
plaquette. In other words, we can write £^ — i^J^^ where 
Jfj, is the vortex current of a 2+1 dimensional XY model. 
Assuming the latter is described by the critical point of 
the 0(2) ip 4 field theory in 2+1 dimensions, we can-com- 
pute the correlations of Jn by Halperin's analysisEl, as 
generalized in Appendix [D| to the critical point. A com- 
putation very similar to that leading to (3.7) now yields 



c v<liV {k) = (Mk)M-k)) 

= (fc 2 ^ - fc M fc„)(& + Z 2 k 1+ ^ + ...), (3.16) 

where now we choose to define (see Appendix |S|) 
dim[J A1 (r)] = (6 + i] v )/2. The transverse nature of the 
right-hand-side is a consequence of the conservation of 
the vortex current, which is now expressed by the local 
law 



9/j,O v>l 



0. 



(3.17) 



If (3.17) was not violated by an anomaly, then we would 
expect that di m] t 2 tt ( r)] = 2 or r) v = —2. Inserting the 
latter value in (|3.16|) we see that the singular contribu- 
tion from the scaling theory dominates the non-singular 
piece proportional to ci . This is the opposite of the situa- 
tion in the SU(2) symmetric case considered in the body 
of this section. In the present case, we do not expect 
an anomaly in the conservation of vortex current as the 
conservation law is dominated by the scaling theory con- 
tributions, and so the result rj v = — 2 should be ex act in 
3. Integrating out the fluctuations using (3.16) 



D 



and this value of r/ v , we conclude that 



exp 



C2 

2 



E ^a A1 (-fc)(fc 2 <5 A1 - k^k v )a u (k)j 

(3.18) 



for t he cas e of Zxjm- The strong contrast between ( 3. IS ) 
and ( 3.15 ) should be noted. The present resul t ( 3. IS ) 
is midway between the a-Meissner result (|3.15|) which 
expels flux, and the "a-Maxwell" result ( p!l|) which 
allows penetration of a M flux. A more careful analysis 



is required to understand the consequences of ( 3.18 ) and 
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this will be taken up in Section [v| We will find there that 
the BC order can either present, absent, or critical at the 
point at which U(l) spin-rotation invariance is restored; 
the situation depends upon the details of the model be- 
ing studied. For our purposes here, it is satisfying to 
note that, unlike the SU(2) case, the present methods do 
not place the U(l) critical point in the a-Meissner phase, 
which would have required BC order. 



2. Remarks 

We conclude this section by a few remarks on SU(2)- 
invariant systems, motivated by the previous subsection 
on U(l) symmetry. It is interesting to observe that if the 
SU(2) symmetric spin system had Ch{k) ~ k 3 (the value 
suggested naively by ( |3.3|) ), then its W{a tl ) would have 
had a powe r of k in the exponent which was the same 
for systems with U(l) symmetry; in 



as that in (3.1 



such a situation one could envisage a direct second order 
phase transition between the HN and BC phases, as had 
been conjectured in earlier work by one of uscj. We have 
argued above that this behavior of W(a ll ) is not present 
at the critical point described by the 0(3) symmetric (p 4 
field theory. This paper does not find a clear-cut scenario 
under which this behavior could arise generically, as we 
will discuss further at the end of Section VC 2. 
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FIG. 4. Phase diagra mji of the S = 1/2, one- dimensional 
antiferromagnet in ( |l.l| ). The BC and I phases are as in 
Table and have a gap to all excitations; these phases have 
gapped excitations with S z = ±1/2 but are not fractionalized 
in the two-dimensional sense — the S z = ±1/2 excitations con- 
fine immediately in pairs when such one-dimensional chains 
are infmitesimally coupled together to form a two-dimensional 
lattice. The Tomonaga-Luttinger (TL) phase is special to one 
dimension: it is gapless and has power-law spin correlations. 
There is SU(2) spin symmetry along the line f = 1. 



A. Results from bosonization 



IV. ANTIFERROMAGNETS IN ONE 
DIMENSION 

There is now a rather complete understanding of S = 
1/2 quantum antiferromagnets, like those in (1.1), in one 
spatial dimension. This is largely due to the powerful 
bosonization technique, which not only allows classifica- 
tion of various gapless critical states, but also describes 
the renormalization group flow between gapped states 
with numerous possible broken discrete symmetries. We 
will use these well-established results as a convenient lab- 
oratory for testing the methods developed here for two- 
dimensional antiferromagnets. As we show below, when 
adapted to one dimension, our methods reproduce all the 
important universal results obtained by the bosonization 
method. This agreement significantly boosts our confi- 
dence that our lattice models and duality methods are 
also successfully capturing the physics of antiferromag- 
nets in two dimensions. 

O ur dis cussion is divided into two subsections. In Sec- 
tion IV A, we recall the phase diagram of a S = 1/2 
quantum antiferromagnet in one dimension obtained by 
the abelian bosonization method. 



we 



In Section [IV B 
show that an essentially identical phase diagram is ob - 
tained by an analysis of the partition function Z in (1.13) 
when applied to one dimension. 



A convenient frustrated S = 1/2 antifecmmagnet in 
one dimension is that studied by Haldaneo with first 
and second neighbor exchange interactions: 



Ji (i 



± J2&j ■ Sj+2 



(4.1) 



We have also introduced a uniaxial anisotropy parameter 
C so that the spin symmetry is SU(2) only at £ = 1, and 
is U(l) otherwise. Both exchange constants are positive, 
Ji,2 > 0, so the second neighbor coupling is frustrat ing. 

A review of the abelian bosonization analysis of (4.1) 
may be found in Rcf. The resulting phase diagram 
is shown in Fig ^. We exclude the region of large J<xj J\ 
where additional phases with a time-reversal symmetry 
breaking, spin chirality order parameter coexisting with 
BC or TL order are foundEl There are two phases with 
a gap to all excitations, the BC and I, and these are as 
described in Table §. The I phase occurs for large £, when 
the spins clearly prefer to be polarized along the ±z direc- 
tions. The BC phase appears for moderate values of the 
frustration, J2/J1, and the spatial pattern of the bond 
order, Qjj+i, is sketched in Fig 0^-.this phase contains 
the special Majumdar-Ghosh pointES £ = 1, J2/J1 = 1/2 
where the fully dimerized wavefunction, consisting sim- 
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ply of products of singlet pairs of nearest neighbor sites, 
is the exact ground state. 

The Tomonaga-Luttinger (TL) phase is gapless, and 
special to one dimension. It is the analog of the XY 
phase in two dimensions, but instead of exhibiting true 
long-range order, spin-correlations in the x-y plane have 
a slow power-law decay, 



(Six ^j-\-n,x ^iySj+n.y) 



(-1)" 

\ n \nxY 



(4.2) 



for large |n|, where the exponent r\xY < 1 varies con- 
tinuously within the TL phase (rjxy = 1/2 at £ = 0, 
J2/J1 = 0) . Similarly, the correlations of the Ising order 
(—iySj )Z decay with an exponent rji = I/t/xy, while 
those of the BC order decay with an exponent also equal 
to 77/. 

At boundary between the TL and BC phases, and also 
at the boundary between the TL and I phases, we have 
rjXY = T)i = 1; actually, there is also a logarithmic cor- 



rection to (4.2) at these boundaries-we do n ot dis play 
this here but note that our methods in Section [VB| also 
reproduce this correction. Within the TL phase we have 
t]xy < 1- The boundary between the BC and I phases is 
also critical, and has t]xy > 1- 

The manner in which the BC and I orders vanish at the 
boundaries of their respective phases is also reviewed in 
Ref [34f both order parameters vanish with an essential 
singularity at their boundaries to the TL phase, while 
they vanish with the same continuously varying power- 
law on opposite sides of the BC-I boundary. 



B. Dual lattice models 

We will now describe the properties of the lattice model 
Z in (1.13) adapted to one dimension: this latter model 
will be denoted We will see that the phase dia- 

gram of a generalized class of such models is identical in 
structure to that of Fig ^. 

First, let us explicitly write down the form of Z^; we 
have 



Z w = £ / II dn 3 S(n 2 3 - 1) exp (-S u - S$> - Si 

(4.3) 

In present situation in two spacetime dimensions, the site 
indices j and j extend over the square and dual-square 
lattices respectively, and the Greek indices /i, v, A . . . 
extend over x, r, and the parameter r)j in the Berry phase 
term is now (— l) Jjl . Furthermore, t he A flux can point 
in only one direction, and so Sa in ( |1.11| ) is replaced by 



S 



(4.4) 



so that the integer qj is simply a scalar and carries no 
spacetime vector index. We wish to map the phase di- 
agram of Z^y in a parameter space where the coupling 



in S n in (1.9) is allowed to acquire a uniaxial anisotropy 
and reduce the spin symmetry to U(l). 

We proceed with the duality transformations as in Sec- 
tion H Then the expression for the partition function in 
(|2.4|) is replaced by 



(4.5) 



where now the integers cij carry no spacetime vector in- 
dex, and the action S a is replaced by 



si" 



3 



a°) 2 , (4.6) 



TjjS/j.r instead of (0 



with a? now satisfying e^a 
a convenient choice has a° taking the values 1 and on 
alternating spatial columns i.e. a? = (1 + (— 1)^)/2. 

We now need to underst and the consequences of the 
integration over the rij in (4.5). Because our purpose is 
to compare with the properties of the model (4.1) which 
has a uniaxial anisotropy and only a U(l) symmetry, we 
choose to introduce a similar anisotropy in S n . Actu- 
ally it is convenient to take the easy-plane limit, and to 
assume that the n fields are completely localized in the 
x-y plane in spin space. In such a limit the integral over 
the n fields in ( |LE| ) can be evaluated completely, and a 
precise mapping to a dual lattice model is obtained. We 
will then examine at the general structure of the resulting 
dual lattice model, and discuss how its parameter space 
allows us to relax the easy-plane restriction, and even 
find regimes in which the spins are localized in the ±z 
directions, and also special parameter values for which 
full SU(2) spin symmetry is restored. So, in the easy- 
plane limit, we parameterize iij in terms of an angular 
field Gj as noted in (1.14). To aid the integration over 



the Qj , we perform the usual innocuous change of writing 
the action S n in the periodic Gaussian (Villain) form 



£ n ,u(i) = iz- £ (A M 0j - 27rm iM ) 2 , 



2.9^ 



(4.7) 



;;■/"< 



where the 



are the integers on the links of the di- 



rect lattice which ensure the periodicity of the action in 
the angular Qj. These integers also have the utility of 
allowing us to measure the vorticity within each plaque- 
tte: this is simply e^^A^m^v- Now recall the discussion 
in Section pn in the paragraph containing (3.16), where 



we noted that in the easy-plane limit the A flux through 
a plaquette is simply 7r times the vorticity. U sing this 
result, and the modified from of the action in (4.7), the 
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partition function Z^ 1 ' transforms to the following form 
in the easy-plane limit 



Z V(i)= E |lI J J ex p(" 5 ^(i) 



s a ] - iTr^e^uajA^mjp) . (4 



Before proceeding to a duality mappin g o f (4.8), we 
present an alternative representation of ( j4.8| ) which ex- 
poses its close relationship to the Zi gauge the ories of 
Refs. H,p7|- We notice that the last term in (4.8) is only 
sensitive to whether the integer rrij^ is even or odd. So 
we write 



= Vn 



1- Sj 



3+H 



(4.9) 



where /j M is an integer vector field on the direct lattice, 
and Sjj+p = ±1 is an Ising variable on the links of the 
direct lattice; notice that Sjj+ P is unoriented (unlike rrij p 
and fjp) and this is reflected in the choice of notation. 
At this stage, it is now easy to perform the summation 
over the Oj independently on every dual lattice site, and 

to obtain the alternative representation of Z 



(i) 

U(l) 



{s jtj+[l =±l}{f jl >} J 3 \ o n 

2^fA^ l d j 7T 

y 3,fi v 



'3,3+H 



4E ,:I *<•<•-• 



(4.10) 



The first term in the action of this form of ^u(i) * s * nc 
standard Maxwell term of a Z^ gauge theory, with Sj,j+/i 
the Zi gauge field. The coupling K is related to e 2 by 



E 



-e 2 n 2 /2 



2if _ n=-oo 



(-l)"e- e2 " 5 



(4.11) 



/2 



n— — oo 



if is a monotonically decreasing function of e 2 , with 
K = 7r 2 /(4e 2 ) - (ln2)/2 as e 2 ^ and K = 2e~< ? I 2 



as e — *■ oo. The last term in ( 4. 1C ) is the Berry phase 



of the Z2 gauge theory: this is identical in form to 
that in Ref. || and was obtained after using the iden- 
tity e^Aydj = r)j5^ T . The connection to the models of 



Ref. |8J becomes clearer if we write ( 4.1C ) using the sim- 
ple cosine interaction of the angle 9j, rather than the 
periodic Gaussian, Villain form: 



J u(i) 



e /n«»i«p(*Eii 



{ S J,3 + A= ±1 } 



^ 



„ ^ Sj,3+H COS 
^ 3,fi 



This is the action of an XY model coupled to a Z 2 gauge 
field, with an additional Berry phase term associated 
with the 5=1/2 nature of the underlying antiferromag- 
net. Note that the angular variable in the XY coupling 
is 0/2, whic h is half the angle determining the orienta- 
tion of n in ( 1.14 ): this half-angle variable is of course the 
reason for the appearance of the Z2 gauge degrees of free- 
dom. Thus we have established the rather surprising re- 
sult that simply by i ntrod ucing an easy-plane anisotcopy 
in the model Z in ( 1.13 ) of Sachdev and Jalaberto we 
obtain the XY-Z2 gauge model of Senthil and FishcrtJ; 
this result holds also in two d imensions, as we will see in 
Section [y|. The model (4.12) generalized to two dimen- 
sions, but without the Berry phase, was studied recently 
by Sedgewick et alSl. 

We now return to the original expression for zjiL in 



(4.8) and proceed with the duality mapping: this will al- 



low evaluation of the 6j integrals and will also write it 
in a form with only positive weights. First, we perform 
the summation over the m JAI by using the following gen- 
eralized Poisson summation identity on every link of the 
square lattice 



exp 



— — (x — 2irm) + imy 



E ex p 

j=— 00 



(4.13) 



where to, J are integers, and C > 0, x, y are arbitrary 
real numbers; this replaces the sum over the m JM by 
a summation over a different set of integer "currents" , 
Jj M , also residing on the links of the square lattice. The 
advantage of this form is that the integrals over the 0j 
can be performed independently, and yield only the con- 
straint that the Jj P currents are divergenceless. In this 



manner, Z^) x in (4J3) is shown to be exactly equivalent 



to 



Z u(i) = E ex P 



Ek-4) 2 



3,3' 



Jj'u + 



, (4.14) 



where the prime on the summation indicates the con- 
straint 



= 0. 



(4.15) 



The mapping from (4.8) to (4.14) yields only a on-site 
coupling between the currents with 



5jj> Sf+v 



(4.16) 
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(this KjZ is not to be confused with the Ising gauge 



coupling K in (4.10, 4~12] ); the latter does not carry any 
indices). However, to reproduce the full phase diagram 
of the antifcrromag net #W in ([□]) it is necessary to 

""""As 



allow for additional near-neighbor couplings in 



will become clear from our discussion in Section [V B 1 



below, these new terms allow us to include the couplings 
between the z components of the rij which we re d ropped 
when we specialized the SU( 2) i nvariant S n in ( |l.9[ ) to the 
easy-plane limit 5 n ,u(i) in (4.7). With the off-site terms 



in K?p we will be able to restore the SU(2) symmetry 
along special parameter lines of Z^\y and even obtain 
regimes where the spin anisotropy is easy-axis, and the 
spin-ordering is Ising-like. 

The expression (4.14) is one of the final dual forms of 



AD 



■^u(i) - man y physical properties will be most transparent 
in this formulation. However, for future technical pur- 
poses, it is useful to perform a few mor e manipulations 
on it. We can solve the constraint ( 4.15 ) by writing 



J, 



(4.17) 



where pj is an integer on the sites of the dual lattice. 
Then, after defining £j = 2pj + etj, it is easy to see that 
the summation over the aj c an b e performed explicitly. 
In this manner we find that (4.14) reduces exactly to 



Z u(i) = ex P 



3 



(4.18) 



where ctj is a fluctuating Ising variable which denotes 



whether £j is even or odd, 



ctj = 1 - 2(£j mod 2), 



(4.19) 



while Ej = (— 1)^ is a fixed oscillating field which arises 
from the background a? field, and is ultimately related 
to the spin Berry phases. The coupling Kd is the "dual" 
of the coupling K 



tanh Kd 



-IK 



(4.20) 



where K was defined as a function of e 2 in (4.11); note 
that Kd is a monotonically increasing function of e 2 , with 
K d = 2e~ w2 /( 2e2 ) as e 2 -> 0, and K d = e 2 /4 - (ln2)/2 
for e 2 — > oo. 



1. Observables and limiting cases 

We pause in our chain of duality mappings to discuss 
the physical interpretation of the new degrees of free- 
dom we have introduced. The physical picture becomes 
clearest after an understanding of the interpretation of 



the current s Jj a in (4.14). Analysis of the duality map- 



ping from Q to (pT|) along the lines of the dualities 
of quantum XY and boson models in Refs. |3^j3^ shows 
immediately that these currents are precisely those of a 
boson representing the spin-fli p op erator. In other words, 
if in the original spin model (hi) we make the identifi- 
cation 



J 3> x 



iS 



j,y 



(4.21) 



and perform the usual boson-duality transformations (as 
discussed by S0rensen et al£3) to a model of interacting 
integer- valued current loops, the will be the spatial 
and temporal (i.e. the density) components of the cur- 
rent of the boson in ( [4.21 ) . More precisely, we can couple 
an external chemical potential to the boson charge and 
thence deduce that 



S 



3,z 



f A /- = A /- 



(4.22) 



So the spatial gradients of the £j are related to the value 
of Sj. z , and a phase with I order will have a staggered, 
zigzag pattern of the £j. 

The above interpretation can be reinforced by exam- 
ining the limiting case e 2 — * oo (or equivalently, Kd —> 
oo). Then in (4.14), the aj fluctuations are quenched to 

(-l)^V- With this value ' 



2 j = a j and so e M „A,. 



the partition function in (4.14) is seen to be precisely t 



boson current loop model examined by Otterlo et al£J 
(but in two spatial dimensions) at a mean boson density 
of 1/2 per site. This last density is, of course, precisely 
that needed to have a state with (J2j Sjz) = in the 
spin model. With the purely on-site coupling in ( 4.16| ) 
this boson-loop model is expected to always be super- 
fluid: in one spatial dimension the superfluid (or XY) 
correlations are only power-law, and so the spin system 
is in the TL phase. 

With t he interpretation of the spin operators in 
( 4.21 . i~2^ ), we have just seen that at Kd — oo each site 
has an average density of 1/2 spin-flip bosons. This is 
just as expected from the usual microscopic interpreta- 
tion of the original spin antiferromagnet. However, what 
about the case when Kd < oo ? In this picture, there now 
appear to be fluctuations in the mean boson site density, 
e TV A v aj/2. Where do the missing bosons go ? The an- 
swer is that they reside on the links of the lattice; indeed 
if the spins on sites j and j + x form a spi n-singlet, this 
is expressed in terms of the bosons in (4.21) as the single 
boson state 



|0) 



(4.23) 



where |0) is the reference state with both spins down. So 
a boson on a link corresponds to a spin-singlet bond, and 
each site on the ends of the link will lose half a boson. So 
our present model for spin systems goes beyond that of 
Otterlo et alEB by allowing easily for coherent occupation 
of bond orbitals. This sets the stage for states with BC 
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oo 



FIG. 5. Phase diagram of the one spatial- dimen sional 



model Zy)^ for the case of the on-site coupling (4.16). The 



J u(i) 



phase boundary is in the Kosterl itz-Thouless universality 
class, as discussed in Section [VB2. 



order, and the present interpr etat ion is seen to be com- 
patible with identification in (2.6) of the "bond-charge" 
density with the energy associated with the aj fluctu- 
ations on the dual lattice site in the center of a bond. 



Using the variables appearing in ( 4.18 ) , it is evident that 
the bond modulation is now measured by 



(4.24) 



where j resides on the middle of the plaquette between 
j and j + x. With this identification, we can see that 
BC order is present in Z 1 ^^ in the limit g —> oo. From 



(4.18), we see that l T is independent of J, and therefore 
so is (Tj. From ( 4.24 ), a constant Cj leads to an oscillating 
Qij with BC order, with the values aj = ±1 identifying 
the two dimcrized states. 

We now have enough information to construct the 

phase diagram of ^u(i) ^ or ^ ne case °^ ^ ne P ure ly on-site 



coupling (4.16). From (4.18) it is clear that there can 
be no phase transition at g — 0. Moreover, at Kd = 0, 
( [4.18 ) reduces to a height model in 1+1 dimensions, and 
this must have a roughening transition at a critical value 
of <?; this will be discussed further in the following sub- 
section. These results then allow us to sketch the phase 
diagram in Fig 0. Note there is no I phase; this will ap- 
pear below when we perform generalized analysis which 
allows for further neighbor couplings in the K^,. 



2. Field theory and phase diagram 



We now perform a field theoretic analysis of ^un.) 



al- 



idea, as mentioned above, is to view the partition func- 
tion as a model of interfacial "heights", £j, in 1+1 di- 
mensions. Then the K^, control the roughness of the 



interface, while the first term in (4.18) indicates that the 
heights are preferentially even (odd) integers on every 
even (odd) spatial column. This height model can be an- 
alyzed by entirely standard technology, which is reviewed 
in Appendix [5|. The final result is that Z^)^ is related 
to the continuum sine-Gordon model in 1+1 dimensions 
with action 



c(l) _ 
D sG 



fi'r 



(|(a M o ) 2 -Fcos(2^o)) , (4.25) 



where r is a spacetime co-ordinate. In terms of the con- 
tinuum field 0o i the I order is measured by sin(7T0o), 
while the BC order is measured by cos(7r</>o)- The cou- 
pling Y is related to the parameters of ^uA) (E5). 



The action ( 1.25 ) is of the standard sine-Gordon type, 
and its properties can be read off from existing results. 
For small g, the coupling Y scales to and its only ef- 
fect is a renormalization of g to g. This yields the TL 
phase, and it is stable provided g < 2ir. It is gapless 
and has power-law decay of correlations, controlled by 
the fluctuations of the free field 4>o- With the identifica- 



tion of the BC and I order parameters below (4.25), we 
can now compute that they both decay with the expo- 
nent r\i = 2n/g. To compute the corr elations in the x-y 
plane in spin space we note from (4.7) that fluctuations 
of the angular variable 9 are controlled by the continuum 
action 



S 



(i) 

XY 



1 

2f) 



d 2 r(d^) 2 ; 



(4.26) 



Now computing the correlat or o f e 1 ® we deduce that x-y 
spin correlations decay as (4.2) with rjxy = fl/(27r) = 
1/rji. Note that all of these exponents are identical to 
those o btaine d by abelian bosonization, as reviewed in 
Section tVA. 



lowing for a general short-range coupling K^., . The main 



The TL phase is unstable for g > 2tt, and the ground 
state acquires a gap induced by the flow of \Y\ to large 
values. For Y > 0, the field 0o is preferentially pinned 
at the values (f>o — 0,1, and these are easily seen to cor- 
respon d to the BC phase; the case with on-site couplings 
as in (4.16) is therefore seen to correspond to Y > 0, 
as (4.25) describes the phases in Fig || Similarly, for 
Y < 0, we have <f>o — 1/2,3/2, the values correspond- 
ing to the two I states. These results imply the phase 
diagram shown in Fig ^. There is complete consistency 
between all the universal properties of Fig ^, and those 
obtained by abelian bosonization for the model of Fig |J. 
Establishing this was the primary purpose of this section. 

A notable feature of Figs [| and || is that SU(2) spin 
symmetry is restored at the boundary between the TL 
and I phases. Moving along this line we see that the 
ground state initially has SU(2) symmetric spin correla- 
tions decaying with an exponent t/xy = Vi = 1- How- 
ever, there is eventually a continuous phase transition to 
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. TL 



FIG. 6. Ph ase d iagram of the one spatial-dimensional 



model S^q in (4.25). The same phase diagram also applies 



to ^u(i) m (4-14) with general off-site couplings beyond the 



on-site coupling in ( 4. j_6 ) . The topology of this diagram is 
identical to that of Fig ^, as are all the universal long-distance 
properties of the TL phase, and of the phase boundaries. 
There is SU(2) spin symmetry along the grey hatched line. 
This line includes the boundary between the TL and I phases, 
but does not represent a phase boundary within the BC phase. 



the phase with BC order. So we have succeeded in un- 
derstanding the phases of the SU(2) model via a detour 
into models with only a U(l) symmetry. 

We conclude this section by noting that closely related 
results are obtained in an analysis based on the CP 1 
formulation presented in Appendix [X]. This extension is 
discussed in Appendix [|J This is in keeping with our 
contention that the n field formulation of Section | and 
the CP 1 field formulation in Appendix [A| have similar 
phase diagrams. 



V. U(l) SYMMETRY IN TWO DIMENSIONS 

Bolstered by the success of our approach in one di- 
mension, we forge ahead to the application of the same 
methods in two dimensions. We are interested in mod- 



els of two-dimensional antiferromagnets like (1.1), but we 
now allow the spin exchange terms to acq uire a uniaxial 
anisotropy like that in the first term in (4.1). We will 
carry out duality mappings on such spin systems with 
U(l) symmetry, and then search for lines in generalized 
parameter space where SU(2) symmetry can be restored. 
This should considerable light on the properties of SU(2)- 
symmetric quantum antiferromagnets. 

O ur ini tial analysis follows the same steps as in Sec- 
tion IV B; the only change is that some of the fields 
acquire additional vector indices, associated with the 
change of spacetime dimension to 3. We in trodu ce spin 
anisotropy into the partition f uncti on Z in ( LIS ) by pa- 
rameterizing the n field as in (1.14) and replacing 5 n in 
(1.9) by iSn,U(i) m (4/7). Then proceeding with the dual- 
ity mapping at the beginning of Section O, we obtain the 
following U(l) symmetric version of the partition func- 
tion Z in (p.4|); this generalizes (4.8) to two dimensions 



Z U(1) = £ /lI d ^' ex p( _ ' S 'n,U(l) 

(n,>,mji.} 3 

- S a - in 2J eftuxa^Aurrijx 



(5.1) 



As in Section IV B , we can obtain alternative represen- 
tations of (5.1) which expose its close connection to a Z2 
gauge theory. Indeed, these representations for Ztj(i) are 

identical to those for zjiL in (|4.10 ) and ( 1.12 ) with the 
only change being that the summation over fj, extends 
over the three spacetime directions: this al terna tive for- 
mulation of ^u(i) was presented earlier in ( |1.15 ). 

Next, we proceed with the two-dimensional generaliza- 
tion of the dual ity mappings discussed below (4.12); we 
generalize (4.14) to 



Z\j(i) = ^2 exp 

{ a jfj. »^j>} 



-t£(%*-<J 2 



"■jj> I "in • r, II J j'y -t- 9 



3,3' 



(5.2) 



where, as before, the prime on the summation indicates 
the constraint (4.15). The particular model (5.1) corre- 



sponds to the purely on-site values of the coupling Klf., as 



in (4.16). However, as discussed in Section IVB, we find 



it useful to consider models with near neighbor ifjtv to 
allow for phases with easy-axis anisotropy and parameter 
values with enhanced SU(2) symmetry. Another useful 
form of ( |5.2| ) is obtained by solving the constraint ( 4.15| ) 
by 



(5.3) 



with pj\ integer; this generalizes (4.17). Wi th th is pa- 
rameterization, the two-dimensional form of (4.18) is 



J TJ(1) 



Yl exp 



3,3' 



3 



3,3+P- a 3,3+fL 



(5.4) 



where = + a jp i s an integer valued vector field on 
the links of the dual cubic lattice 
degree of freedom which is pinned to the ij^ by 



and <7j T-j-n is an Ising 



J JJ+A = 1 _ 2 ( 



mod 2), 



(5.5) 



generalizing (4.19). Notice that, unlike Ij^, o~j,j+ti is un- 



oriented, and has only a single value ±1 on each link; this 
is the reason for the different choice of notation for these 
two fields. The fixed field is likewise unoricntcd 

and takes values ±1. It arises, as in (4.18), from the 
Berry phase terms in •Z'tj(i) > anc i its values must satisfy 
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n 

spatial □ 



3,3+f* 



(5.6) 



about every spatial plaquette. All other plaquettes have 
Y[ D = 1< The coupli ng K d is a mono tonic function 

of e 2 given, as before, by (4.20) and (4.11). 



A. Observables and limiting cases 



Continuing our parallel to the discussion of one- 

we gener- 



dimensional antiferromagnets in Section IV B 



alize the discussion of Section IV B 1 to two dimensions 



The generalization of the observables is immediate. The 
spin-flip site boson density is now measured by the I flux 
piercing the associated plaquette on the dual lattice: 



S 



3^ 



(5.7) 



(generalizing ( 4.22 )). The BC order in ( 4.24 ) needs some 
notational modification: 



13, 3+% ~ £ 3,3+y <J 3,3+Vi 



(5.8) 



where the direct lattice link on the left-hand-side inter- 
sects (after projecting the r co-ordinate by half a lattice 
spacing) the dual lattice link on the right-hand-side. A 
similar relation holds with x <-> y. 

Next, we turning describing the physics of the limiting 
parameter values in the phase diagram of Zyny We will 
see that there is much more structure than that in the 



corresponding discussion in Section IV B 1 



As before, we first consider the limit Kd 



Now 



the results are essentially identical to those in one di- 
mension: the values of asj^ are pinned to a?, and so 
e^xA^a^x = Vj^^T- The partition function ( |5.2] ) is ex- 
actly the boson current-loop model studied by Ottcrlo et 
al£3. With the purely on-site coupling K^, in (4.16), 



they found that the boson-loop model was always a su- 
pcrfluid, which corresponds to a spin system with long- 
range XY order. With a nearest-neighbor coupling in the 
Kjj 1 /, Otterlo et al. found states with a boson site density 
wave, which corresponds here to spin states with I order. 
In the limit g — > oo, Zym reduces to a non-trivial 



model, unlike the case in Section IVB1. From (5.4) we 
deduce that the integer vector field £j M must be curl- 
free, implying it can be expressed as the gradient of an 
integer-valued field on the sites of the dual lattice. The 
relationship (5.5) then implies that Ising field Oj,j+/; is 
also 'pure gauge' i.e. 



a 3,3+ij- ~ ^3^3+V 



(5.9) 



where <rj = ±1 is an Ising field on the sites of the dual 



lattice. So 



J U(1 



(9 




(5.10) 



This is the partition function of an ordinary Ising model 
on a cubic lattice, with a nearest neighbor exchange of 
modulus Kd, and the signs of the exchange chosen so 
that every spatial plaquette is frustrated. The expres- 
sion (5.10) could also have been obtained by appl ying the 
standard Ising duality to the Z^ gau ge the ory in ( LIE ) in 
the limit g = o o; the Berry phase in ( 1.1 5| ) appears as the 
frustration in ( 5.10 ). Als o no te that in the framew ork o f 
the parameterization in (5.2), the steps leading to 



correspond precisely to the discussion in Section |I that 
the large g limit in t he pr esenc e o f U(l) symmetry im- 
posed the constraint (2.10) on ([2~g|). The Ising model in 
fl5.10| ) was studied in Ref. 16|,[l8] with the motivation also 
provided by S = 1/2 quantum antiferromagnets. The 
context of U(l)-symmetric models is quite different here, 
and the earlier motivation is close to the SU(2) symmetric 
models to be considered in Section VI. From this earlier 



work on this frustrated Ising mod el, a nd from methods 
to be described in detail in Section VB , wc conclude that 
the "ferromagnetic" phase found at large Kd has BC or- 
der. Also, the Ising "paramagnetic" phase at small Kd is 
fractionalized (F) with no broken lattice symmetry. This 
F phase had no analog in one dimension, where the large 
g limit yielded only the BC phase. 
Turning to the limit Kd 



0, ^u(i) m (5-4) becomes 



a model of interacting /iaZ/-integer currents with short- 
range interactions. After a rescaling of currents and cou- 
pling constants, this model is formally equivalently to one 
with integer currents, which correspond in turn to boson 
models at integer filling studied in Rcfs. ||^,|39|. These 
have a transition between a superfluid phase at small g 
and an insulating phase at large g which is in the univer- 
sality class of the 0(2) (p 4 field theory in 3 dimensions. 
For the spin model, the superfluid phase corresponds to 
the XY phase, while th e insula ting phase is F phase. 

Finally, as in Section IV B 1 , the g ^ limit is trivial 
and the ground state always has XY order. 

Putting all these results together, we can sketch the 
phase diagram in Fig [?] for ^u(i)- Actually, determin- 
ing the structure of the phase boundaries in this figure 
required some additional field-theoretic analysis which 



appears in Section VB 



B. Field theories 



We now proceed to a field-theoretic analysis of the 



model -Zu(i) in (5.4) following the same method used to 
^ Qo„f;™ It^roI T he technical details of 



analyze Z^^ in Section PVB 2 



the derivation are similar to those used in the one di- 
mensional case in Appendix [|| and their generalization 
to two dimensions is described in Appendix [g|. 

The field theory obtained in this manner turns out to 
be very closely related to that, studied recently by Lan- 
nert, Fisher and Senthil (LFSp, although in a somewhat 
different physical context. They were studying charge 
fluctuations in a correlated superconductor using models 
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FIG. 7 . Ph ase dia gra m of two spatial-dimension model 



HJ(l) 



in (1.15) or (5.2) or (5.4) for the case of the on-site cou- 



pling (4.16). The thin lines are second-order phase boundaries 
while the thick line is a first-order boundary. The field theo- 
ries for the p hases and the critical properties are discussed in 
Section |VB. 



of in tera cting boson Cooper p airs 
tion 



As we noted in Sec- 



[A 



and again in Section IV B 1 , our antiferromag- 
netic spin models can also be interpreted as interacting 
boson Hamiltonians, but the boson is now a spin-flip op- 
erator, as in (4.21). Casting aside the different physical 
interpretation of the bosons, the results of LFS can be 
directly applied to our models: we simply have to map 
their superfluid state to a XY phase of antiferromagnets, 
their boson site charge-density-wave to an I phase, while 
the interpretation of BC order remains the same in both 
cases. We will follow the same general approach as LFS, 
but will extend their results in several directions to ad- 
dress the questions of interest here. 

In one dimension, as presented in Section IV B 2 , the 
field theory describing the phases of the antiferromagnet 
was the sine-Gordon model in ( 4.25| ) . In two dimensions, 
the analogous procedure in Appendix ^| shows that the 
central actors are poi nt vo rtices in the wavefunction of 
the spin- flip bosons in ( 4.2l| ) . We will denote the complex 
field which creates or annihilates single vortices by 0(r), 
where r is a spacetime co-ordinate (the boson wavefunc- 
tion winds by a phase of 2tt around such a vortex) . Also 
essential for a description of the phases, and emerging 
naturally in the derivation in Appendix |Gj, are the dou- 
ble vortices <f>(r) (the boson wavefunction winds now by 
4tt). As discussed by LFS, and in Appendix ^, the Berry 
phase terms induce frustration in the hopping Hamilto- 
nian for the vortices: consequently, the low energy vor- 
tices have momenta near both (0, 0) and (jr, 0) and it is 
necessary to include both sets in the field theory. We 
denote a convenient linear combination of fields repre- 
senting vortices near these minima by 0i(r) and 4>2{f) 
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Description 



Complex field operators for ±2-7r v ortices in the 
spin-flip boson fa 1 " = S+ (see (fym)). The Berry 
phase frustrates the hopping of such vortices, and 
their dispersion (in a convenient gauge) has min- 
ima near momenta (0,0) and (tt, 0). 0i and <f>2 are 
convenient linear combination of continuum fields 
representing fluctuations near these minima (see 
Appendix O). By the boson- vortex duality, ±27r 
vortices in 01,2 change b boson number by ±1, and 
therefore carry S z = ±1: such vortices must be 
present in both <j>i and (f>2 together to ensure a 
finite energy cost. 



Complex field operator for ±4-7r vortices in the 
spin-flip boson w = S+ (see ( 1.21 )). Its hopping is 
not frustrated and the dispersion has a minimum 
near (0, 0). By the boson- vortex duality, ±2ir vor- 
tices in $ change b boson number by ±1/2, and 
therefore carry S z = ±1/2 i.e. they are spinons. 
The relationship to 0i,2 is $ ~ 4>i4>2 (~ indicates 
that the two sides have correlators with identical 
physical properties). 



Real order parameter for I order, ^i 



Complex order parameter for BC order, Sfsc 



TABLE II. Summary of all the continuum fields intro- 
duced in Section |v] and their physical properties. 



(see Appendix |c|). The motion of the double vortex, 
<f>, is unfrustrated, and it only has a single minimum at 
momentum (0,0). We list these fields, and their physical 
properties, in Table 0; this table also contains additional 
fields to be introduced later in Section V C . 



An explicit derivation of the effective action for 0i, 
02, $ is discussed in Appendix ^. Here we will begin 
by noting the transformations of these fields under the 
symmetries of the lattice Hamiltonian, and then use these 
to eventually write down the most general effective action 
consistent with these symmetries. First, as is always the 
case in the boson- vortex duality mapping, the fields 0i, 
02 5 $ all carry charges with respect to an internal, non- 
compact, U(l) gauge field 6 M (this U(l) gauge field is 
entirely distinct from the compact U(l) gauge fields A^ 
(Section |) and (Section |v| and Appendix [A])). The 
fields 0i,2 carry 6-charge 1, while $ carries 6-charge 2. 
Lattice symmetries lead to non-trivial mappings among 
these fields, as discussed by LFS: 



7r/2 rotation about a direct lattice site: 



e i7r / 4 0i, 



e-™/ 4 02 



x translation by a lattice spacing: 
y translation by a lattice spacing: 



92, 
i02, 
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Label 



BC 
I 

XY 
F 



Characterization 



(I0i| 2 

<01> = ( 



102) 7^0 
-|02| 2 >/O 
02> = <$> = 

W>i) = <^2) = o 



TABLE III. Characterization of some of the state proper- 
ties in Table | in terms of expectation values of the fields (j>i, 
02, and $. States which satisfy more than one property are 
also possible, and obey the union of the corresponding expec- 
tation values e.g. the BC+F state has (0102) 7^ 0, {$) 7^ 0, 
and (0i> = {4,2) = 0. 



(5.11) 

<£> remains invariant under all of the above lattice trans- 
formations. 

Before writing down the effective action, we note the 
connection of these fields to some of the orders charac- 
terizing the antiferromagnet as presented in Table The 
results follow from simple symmetry considerations and 



are summarized in Table [I] 



The BC order breaks only a lattice symmetry, and 
must clearly be neutral under the 6-charge. Examina- 
tion of (5.11) and the symmetries of the states in Fig |l| 
then easily shows that the BC order parameter is the 
composite field 0*02- The same conn ection can also be 
established directly by applying (5.8) to the analysis of 
Appendix The phase of the mean value of this or- 
der parameter distinguishes between the various states in 
Fig l|: arg(0*02) = 0, tt/2, tt, 3n/2 yields the four states 
like Fig pi, arg(0^ 2 ) = tt/4, 37r/4, 5ir/4, 7ir/4 yields the 
four states like Fig |l]b, and all other values correspond 
to an 8-fold degenerate state with co-existing dimer and 
plaquette order. 



For the I order, we note from ( 5.7 ) that Ising ordering is 
associated with a staggered pattern of boson site density; 
in the boson-vortex duality this corresponds to 6-charge 
neutral, staggered, vortex current loops. The structure 
of the vortex eigenmodes in Appendix |g] then shows that 



the I order is |0i 



102 



The XY order corresponds to boson superfluidity. This 
is only possible if none of the vortex fields have con- 
densed, and hence the vanishing of all expectation values 



in Table EI] 



Finally, as discussed by LFS, the F states appear when 
$ is preferentially condensed over 01,2- In the dual boson 
language this means that excitations with boson charge 
1/2 are liberated (these latter bosons are the mselv es ±27t 
vortices in the double- vortex field <&); from ( 4.21 ) we see 
that these are simply the S = 1/2 spinons of a F state. 
Also, it is evident from Table III that co-existence be- 
tween the F and XY orders is not possible, while BC+F 
and I+F states can- in principle, exist. In the approach of 
Senthil and FisherS XY+F states did emerge in a larger 
space of models and degrees of freedom than those con- 



sidered here. Here, we will show that XY+F states do not 
appear in the models considered so far, but are present 
in the models with non-collinear spin correlations to be 
considered in Section fv|. 

The derivation in Appendix |g|, or application of the 
symmetries in ( 5 . 1 1|) shows that action describing the 
phases of ^u(i) is the spacetime integral of the La- 
grangian 

£ = - ibM 2 + \(d^- ibM 2 + - 2^)$| 2 



+ |(e Mi ,A« 2 + S i(|0i| 2 + |0 2 | 2 ) + 
■., -r C.C. + — 1 (|0i| 2 + |02 



^$0*02 



Ml 

2 



2 ) 2 + y|1>| 4 



+ «l|0l| 2 |0 2 | 2 + « 2 |$| 2 (|0l| 2 + |0 2 | 2 ) + M0102) 4 + CC. 

(5.12) 

We have scaled all the field variables to make the co- 
efficients of the gradient terms equal to unity; we have 
also set all velocities equal to unity, and, for simplicity 
of notation, ignored their possibly different values. The 
parameters s± and S2 tune the system across the various 
phase transitions, and we will shortly discuss the phase 
diagram as a function of si and S2- The cubic coupling w 
(no t con sidered by LFS) is allowed by all the symmetries 
in ( 5.11 ), and by charge neutrality: it will play a central 
role in determining important features of the phase dia- 
gram. The various quartic terms couple invariant scalar 
densities and usually only play a peripheral role in sta- 
bilizing the action. Finally, the vs term serves only to 
choose between the different types of BC order. 

The mean-field phase diagram of C has a somewhat dif- 
ferent structure depending upon the sign of v\ . We sketch 
typical phase diagrams as functions of s\ t 2 for these two 
cases in Figs |8| (v\ < 0) and |9| (vi > 0). The sim- 
pler case in Fig || has only three phases: XY, BC, and 
F, and the saddle point always has |0i| = |02| (so there 
is no I order). Any two of the phases can be separated 
by second-order phase boundary (whose field theory and 
order parameter can be easily determined from Table III 
and C). However the phase boundaries necessarily be- 
come first order near the point where the three phases 
meet: this is a consequence of the cubic w term in C, 
as shown recently in Appendix A of Ref. [lO] in a differ- 
ent context but with a field theory with a very similar 
structure. The results of Fig Bwere used to deduce the 
phase diagram of ^u(i) m Fig|7[ Note that the two di- 
agrams have an identical topology. The absence of an 



I phase shows that the purely on-site coupling (4.16) in 



j u(i; 



corresponds to a field theory C with v\ < 0. 



The phase diagram in Fig with V\ > 0, does allow for 
an I phase. So it can only correspond to -Ztj(i) with non- 
on-site couplings Kjj,. A crucial feature of this phase 
diagram is that all phases with I order are well separated 
from the F phase: there is no second order transition be- 
tween the F phase and the I phase (although a strong first 
order transition can always be induced by choosing cou- 
plings judiciously). This feature will play a central role 
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BC 




FIG. 8. Mean field phase diagram of C in (5.15) as a func- 
tion of si and S2 for typical parameter values and Vi < 0. The 
non-compact U(l) gauge field 6 M is ignored in this analysis, 
and the action is simply minimized w.r.t. to mean values of 
the fields </>i,2, As in Fig 0, first-order boundaries are de- 
noted by thick lines, while second-order boundaries are thin 
lines. The BC phase has {4>i,2) 7^ and = \{4>2)\, and 

these can be deduced from the characterizations in Table III 



See also Table IV 



in our dis cussio n of the restoration of SU(2) symmetry 
in Section VC For now, we present a simple argument 
that the separation of phases with F and I order also 
holds beyond the mean-field analysis of C. Let us be- 
gin in the fractionalized F phase. From Table 111, we 
see that the $ field has condensed. The &-gaugc fluctua- 
tions are quenched in such a region (by the usual Higgs 
mechanism), and can therefore be neglected. The onset 
of the I phase involves a confinement transition and so 
one of the fields <pi, <p2 must condense. Let us examine 
the structure of the effective action for ^2 in a region 
where we can safely set ~ and (<&) 7^ 0; from £ the 
quadratic form controlling ^12 fluctuations is 



|9^2| 2 + ( 5 l+«2|(<i>)| 2 )(|0l| 2 + |02| 2 ) 

- - ifl($*)^i0a (5T3) 



Near a second order transition where 0i 2 condense, we 
must diagonalize the quadratic form (5.13) and condense 
the field combinations in the eigenmode with the lowest 
eigenvector. A simple analysis of ( [5. 13 ) shows that, for 
arbitrary ($), this eigenmode has \4>i \ = \<fe\. From Ta- 
ble III, it follows immediately that the condensed phase 



cannot have I order. Therefore the F phase is surrounded 
entirely by with the BC phase (corresponding to conden- 
sation of 01,2) or by the XY phase (corresponding to 
removal of the $ condensate, as is presented in Fig O. 



C. Generalized phase diagrams and restoration of 
SU(2) symmetry 




FIG. 9. As in Fig | but for v 2 > 0. The I phase has 
(4>i ) t^JO, and {4>2) — or vice- versa, as can be deduced from 
Table g. The BC+I phase has (0 ll2 ) / and |(0 X )| / \{<f> 2 }\, 
as indicated in Table tvj. 



The ultimate objective of our analysis of antiferromag- 
nets with U(l) symmetry is to find parameter regimes 
where full SU(2) symmetry is restored, thus enabling a 
description of the phases and phase transitions of SU(2)- 
invariant antiferromagnets. At the corresponding stage 
in our analysis of U (l) antiferromagnets in one dimen- 
sion in Section IV B| we had achieved this objective. We 
found there that boundary between the TL and I phases 
was a critical line, possessing SU(2) invariant correla- 
tions at long distances with a power-law decay (as in 
(4.2) characterized by exponents r\xy — Vi = 1- This 
is, of course, the well-known property of the S = 1/2 
antiferromagnetic Heisenberg chain. Furthermore, as in- 
dicated in Fig ||, the SU(2)-invariant line extends into the 
gapped BC phase, and we obtained the correct and com- 
plete theory of the quantum phase transition from the 
TL to BC phase for SU(2)-invariant antiferromagnets. 

Let us now apply the same strategy in two dimensions. 
As a first step, we have to obtain a parameter regime 
where the ground state has HN order. Unlike the SU(2) 
invariant case in one dimension, the HN state is not crit- 
ical: it has true long-range magnetic order, and an exci- 
tation spectrum of two, gapless, linearly dispersing spin 
waves. By analogy with one dimension, a natural candi- 
date for HN order is the boundary between the XY and 
I phases in Fig 0. However, further consideration shows 
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State 


Condensed fields 


BC 


01, 02, *sc, $ 


I 


0i or 02, */ 


BC+I 


01, 02, *BC, *I, * 


XY 


none 


XY+BC 


*BC 


XY+I 


*/ 


yv l i i ) i 


^ BC i ^= I 


F 


$ 


F+BC 




F+I 


<&, *I 


F+BC+I 





TABLE IV. Exhaustive specification of the states expected 
in the generalized parameter space of C' . Notice that the left 
column contains state labels, and not the label of a property 



III 



that may be obeyed by a state, which was the case in Table 
The right-column contains all the fields that have a non-zero 
expectation value in the corresponding state; if a field is not 
mentioned, it has vanishing expectation value in that state. 
The table above can be deduced directly from the character- 
izations in Table III, and from the field correspondences in 
Table 0. 



that this line clearly does not satisfy the necessary cri- 
teria. This is a critical line describing a second order 
transition: the critical degrees of freedom are the fields 
0i,2 coupled to the non-co mpact U(l) gauge field, and 
the critical field theory is £ in ( 5.12 ) after we set the mas- 
sive field <& = 0. It is expected that all correlators will 
have an anomalous power-law decay along this line, and 
these characteristics are far removed from those of the 
HN phase. 

We have to enlarge the parameter space of our field 
theory to find the HN phase. To this end, note that 



the I (and BC) order in Table [II, were associated with 
composite fields built out of 0i,2- It would clearly pay 
to elevate the I order parameter to an elementary de- 
gree of freedom, as that would then permit co-existence 
of non-critical XY and I order over a wide parameter 
regime, allowing for possible subspace with HN order. 
We therefore introduce the elementary Ising order pa- 
rameter field and, for completeness, the elementary 
BC order parameter field, *£bc- We have now introduced 
all the continuum fields listed in Table [fl], which contains 
a summary of their physical interpretation. 

The relationships of *&/ and ^ bc to the fields 0i.2 
noted in Table || follow immediately from the state char- 



acterizations in Table III . We can now make a catalog of 
possible phases by considering the physical properties of 
all the different combinations of condensation among the 
fields 01,2, ^bc and Wj. This leads to the exhaustive 
listing provided in Table IV. Note the appearance of the 
AF+I phase-this will be our candidate for possessing HN 
order along a subspace in its interior. 

An understanding of the topology of the phase diagram 




XY 




-v"""''iiiiiiiiniimii 

XY+I 



FIG . 10. Schematic mean field phase diagram of £' in 
( |5.14| ) as a function of s± and si with S2 > 0, sbc > 0, and 
vi < 0. The conventions for the phase boundaries are as in 
Fig 0. The grey hatched line represents a possible line along 
which SU(2) symmetry is restored: this line begins within 
the XY+I phase, and then continues along the first-order 
boundary between the XY and I phases; as argued in the 
text it cannot reach the end-point of this first-order bound- 
ary and veers off into another direction not included within 
this cross-section of the phase diagram. 



and the nature of the quantum phase transitions requires 
an effective Lagrangian for all the fields. We ca n use the 
same symmetry considerations discussed near (5.11) to 
extend the Lagrangian C in (5.12) to include \& bc and 
^j. Such a procedure yields 



1 



C = C + -(d^ I f + \d^Bc\ 



si 



4- 



SBC\^BC\ 



WI^I (|0i| 2 - |0 2 | 2 ) - WBC*SC0102 + CC. 



V^BC 



- c.c. + 



(5.14) 



The ellipses represents numerous quartic terms involving 
couplings between the modulus squared of all the fields 
like |*bc| 2 |$| 2 , |*/| 2 (|0i| 2 + |0 2 | 2 ), |*sc| 4 ...; these 
serve mainly to stabilize the action and we will refrain 
from writing them out explicitly. We now have two addi- 
tional parameters, sbc and s/, which can tune the sys- 
tem between its phases by controlling the condensation 
of *5> bc and <J / respectively. The cubic couplings wi and 
wbc follow directly from Table [n], and play a central role 
in determining the phase diagram. 

We can now envisage a four-dimensional phase diagram 
as a function of the parameters s\ 1 s 2 , s/ and sbc- This 
is a space of immense complexity and we will explore 
its structure by describing a few judiciously chosen two- 
dimensional cross-sections. 

Our primary interest is in studying the vicinity of the 
XY+I phase. So we describe the cross-section in the 
si-si plane with s 2 > and sbc > 0- The mean-field 
phase diagram for this case is shown in Fig M. An XY+I 



22 



phase appears, as expected, and we will now study its 
properties in greater detail. 



1. HN order in the XY+I phase 

The XY+I phase only has non-zero, as was indi- 
cated in Table A simple stability analysis about such 
a saddle-point shows that fluctuations of ty'j = — 
$bc, (f>i,2, and $ are all gapped. However, the trans- 
verse component of the non-compact U(l) gauge field & M 
(which is not Higgsed by the neutral $j) yields a sin- 
gle gapless, linearly-dispersing mode: this corresponds to 
the familiar spin-wave mode of the XY order. To obtain 
a state with HN order, we need two spin-wave modes, 
and we can reasonably hope this may happen for special 
values of the parameters. 

As a prelude to demonstrating this more explicitly, let 
us write down the form of the effective Lagrangian that 
will describe the HN phase. In its most symmetric form, 
this is simply the 0(3) </? 4 field theory: 



£* = ~((^«) 2 + ^) + ^) 2 , 



(5.15) 



where a = x,y, z and sh < 0. In the ordered phase, we 
parameterize the fluctuations of ip a by 
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a = {V2Re(y XY ), \/2Im(* X y), (*,) + tfj), (5.16) 



where here (tyj) = \/—6sh/uh, and &xy and ^'j are 
fluctu ating c omp lex and real fields respectively. Inserting 
( 5.16 ) into ( 5.15 ) we obtain the Lagrangian controlling 
the fluctuations about the ordered state: 



C H = \d^xv\ 



\s H \m) 



S H \U H 



6 



4! 



I \2 



21* 



XY 



(5.17) 



Two properties of this form of Ch are especially note- 
worthy. First, there is no mass term for the complex 
field ^ xYi and so there are two gapless spin- wave modes 
at lowest order. Second, a number of cubic and quartic 
couplings are present, but their co-efficients have specific 
algebraic relations betwe en them, tying them to the gap 
in the mode ( y / 2|sjj|) and the single quartic coupling 
Uh'- these relations ensure that the gapless spin wave 
modes survive at all orders in perturbation theory. 



We now search for a n effective Lagrangian like (5.17) 
in C defined in ( |5.14 ). Without loss of gen eralit y, we 
take > and wi > 0. Then, notice from ( 5.14 ) that 
the combination of the wi cubic term and the presence 
of the condensate causes the mass of the <j>2 field to 
be significantly larger than the <pi field. So we proceed 
to integrate out <f>2, along with the double- vortex field 
$ which remains strongly gapped everywhere in Fig ml. 



The resulting effective Lagrangian for <fri and *B'j has the 
4>i controlled by a model of scalar electrodynamics: 

n' ii 1 
C" = \{d„ - tbM 2 + + f L{e^ x d v b x f + S [\^\ 2 + ^|0!| 4 



+^(«) 2 +4(*/) 2 +---> 

where, at lowest order in the non-linearities, 
si =8i- 2wi($i) - 2wi#'j. 



(5.18) 



(5.19) 



We are now ready to make our key nopjtrivial observa- 
tion. We use the well-established resulted that the scalar 
electrodynamics model is dual, in the continuum, to an 
0(2) <p 4 field theory in three spacetime dimensions. Per- 
forming this duality on the (f>i , 6 M degrees of freedom in 
Q5.18 ) , we conclude that C" is equivalent to the following 
dual continuum Lagrangian 



C D = \9^xy\ 2 + s X y\^xy\ 
+ ^(W) 2 + 4«) 2 + 



(5.20) 



where the two-component field of the 0(2) ip 4 field theory 
has been written as the complex field ^xy\ this field has 
a complicated non-local relationship to <fii and 6 M . Fur- 
ther, the parameters sxy and uxy are determined by the 
couplings s' l7 u \ , an d g' in (5.18). The linear dependence 
of s[ on ty'j in ( |5.19 ) implies that the coupling in Zu also 
have linear dependence on ty' T . W ith this dependence, 



the similarity between the (5.20) and 0(3) symmetric 
model Ch in (5.17) becomes evident: the $j d epend ence 
of s'j will lead to cubic terms like those in (5.17). A 
suitable choice of the couplings in £' can clearly ensure 
that its dualized version obeys the same constraints on 
its couplings as those required by 0(3) symmetry in Ch- 
In particular, the mass sxy in Co remains pinned to 
along this 0(3)-symmetric subspace, as is needed to ob- 
tain two linearly dispersing spin-wave modes in the HN 
state. This establishes our claim that HN order appear 
in a subspace with the XY+I phase. 

We have sketched a possible line with SU(2) symme- 
try restored in Fig |l^. A portion of it is in a phase with 
co-existing XY+I order, but it also continues along the 
first-order boundary between the XY and I phases; all 
our above arguments also apply along such a boundary. 
However, within the phases found in Fig there is no 
simple and natural way to extend this SU(2)-symmetric 
line into a non-magnetic phase. If we move above along 
the vertical line si = in Fig 10, we eventually reach the 



second-order boundary between the XY and I phases: we 
argued earlier that this line cannot be SU(2) symmet- 
ric. It is possible to choose parameters such the end- 
points of the two first-order lines in Fig |l^ co-incide: 
then the SU(2) symmetric line can be extended into the 
non-magnetic BC phase, and we have the possibility of a 
second-order transition between the HN and BC phases 
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FIG , 11. Schematic mean field phase diagram of £' in 



(5.14) as a function of sbc and si with si > and S2 > 0. 
The conventions for the phase boundaries are as in Fig (t| 
The grey hatched line represents a possible line along which 
SU(2) symmetry is restored: this line begins within the XY+I 
phase, continues into the XY+BC+I phase, and finally into 
the non-magnetic BC phase; note that sequence of transi- 
tions along this line (moving right to left) is identical to that 
in Fig B (moving left to right), and also to that discussed in 
Section]! 1 1[ Upon going beyond mean-field theory, it is ex- 
pected that tetracritical point M will turn into two bi-critical 
points separated by a short first order boundary between the 
XY+BC and BC+I phases; the line of SU(2) symmetry will 
overlap this first-order line. 



in a system with SU(2) symmetry. But this case is non- 
generic as there does not appear to be any fundamental 
reason for the two end-points to merge. This conclusion 
is also entirely consistent with our earlier considerations 



in Section III 



To find a more natural scheme for a transition to a non- 
magnetic state with SU(2) symmetry we have to search 
the phase diagram in a regime beyond the cross-section 
in Fig [HI This we will do in the following subsections. 



2. SU(2) symmetry with bond order 

A very useful cross-section of the phase diagram is that 
controlled by sj and sbc- This is sketched in Fig [111 in 
the regime s\ > and S2 > 0. As in Section |VC| , the 
XY+I phase is a reasonable candidate for SU(2) sym- 
metry and HN order. Here we discuss the question of 
whether it is possible to continue the line of SU(2) sym- 
metry into the XY+BC+I phase, and then further into 
the non-magnetic BC phase. 

We begin our discussion by the spectrum of excitations 
in the BC phase. This is a confining phase with the fields 
0i,2 condensed (see Table [V). As noted in Table |J, the 



excitations with non-zero total S z are gapped ±27r vor- 
tices in ^1,2 (the relative phases of 0i and 02 arc pinned 
to each other): these vortices have total S z = ±1. In a 
BC phase with SU(2) symmetry this doublet of excita- 
tions must be part of a S — 1 triplet — this is the familiar 
5 = 1 gapped spin exciton expected above a dimerized 
insulator. We will now identify the missing S z = ex- 
citation which completes this degenerate triplet along a 
special line within the BC phase. We will also show that 
condensation of this triplet leads to the BC+HN phase, 
and the latter phase naturally lies along a line within the 
XY+BC+I phase. Thus we will be able to justify the 
entire trajectory of the SU(2) symmetric line in Fig [LT[ 
Within the BC or BC+XY+I phases, the field ^ B c is 
condensed. So we can safely replace ^ bc by (^sc); fluc- 
tuations of ^ bc about this value are massive and do not 
play an important role in our considerations here. Let 
us examine the effective action for 0i,2 and ^ i fluctua- 
tions in phases with (^bc) 0- Cons ider the quadratic 
form controlling 01,2; from (5.12) and (5.14), this has the 
form: 



|d M 0l| 2 + |d M 02| 2 + Sl(|0l| 2 + |02| 2 ) 

- W BC ((*BC>0102 ~ (*BC)^^2 



). (5.21) 



Without loss of generality, we can assume that wbc > 
and (f bc) is real. With these values, the quadratic form 
( |5.21 ) is diagonalized by the linear combinations = 
(0i + 02) /V% and 0' = (0i — 02)/v / 2. The eigenvalue of 
0' is larger than that of 0, and so we integrate out 0' and 
express the effective action £ in terms of and Wj. The 
result has a very simple structure. In particular, there is 
no cubic residual cubic coup ling b etween Vf 1 / and which 
arises from the wj term in (5.14); in terms of 0, 0', the 
term proportional to wj is 4'/(00*' + 0*0'), and upon 



integrating out 0' this leads only to the quartic coupling 
|0| 2 |\l/j| 2 between the modulus-squared of the two fields . 



With these considerations, the new form of £' in ( 5.14 ) 
in the phases with \1/ bc condensed is 



C 



BC 



- *M^| 2 + +^{^vxd„b x f + *|0| 2 + ^|0| 4 



si 



2 

Hi 

2 



;( W + v«) 2 + ^l*'| 4 + ^H*/| 2 - (5-22) 



The S z = ±1 excitations are the vortices in 0, and so we 
should naturally proceed to a duality transformation on 
the 0, degrees of freedom. These fields are controlled 
by scalar electrodynamics and so the transformation is 
similar to that discussed below ( 5.19 ). Scalar electro- 
dynamics is dual to the 0(2) tp 4 field theoryci, and, as 
before, we express this in terms of the complex field $xy 
(which is nothing but the field operator for ±27t vortices 
in 0). So Cbc turns out to be dual to 



£bc = \9^xy\ 2 + +sxy\^xy\ 



^(«) 2 + s 4mr + ^ + v'a^xyi 2 * 2 !. (5.23) 



Ul 
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Now the connection between the £g C and the 0(3) ip 4 
fields theory is plainly vi sible: for sxy — sj and u xy = 
v'j = Auj the theory in ( |5.23| ) is precisely Ch in ( 5.15| ) 
after we identify ip a = (V2Rc(*A'y), \/2Im(* X y), <fj). 
So we have shown that SU(2) symmetry is restored with 
these restrictions on the couplings. For sr = si = sxy > 
0, the theory Ch is in the BC phase and (p a is the field 
operator for the required gapped S — 1 spin exciton; the 
gapped, S z = particle associated with \J// has combined 
with the S z = ±1 vortices in <p to yield the required triply 
degenerate particle. For sh < 0, we have long-range HN 
order induced by the condensation of the S = 1 exci- 
ton; this corresponds to moving along the hatched line in 
Fig [l| from the BC phase to the XY+BC+I phase. The 
critical properties of the transition from the BC phase to 
the BC+HN phase are described by the familiar 0(3) 
field theory in ( 5.15| ). 

The above considerations clearly do not apply to the 
transition from the BC+XY+I phase to the XY+I phase: 
long-range BC order was assumed to be present above, 
and this becomes critical here. The line with SU(2) sym- 
metry plays no special role at this transition, and it de- 
scribed simply by a field theory for ^ bc alone, which is 
obtained by setting all other fields in £ to zero. The 
coupling to the gapless, spin-wave modes of the XY (or 
HN) order is irrelevant. 

All the theoretical justification for the trajectory of the 
SU(2) symmetric line in Fig |ll| has now been provided. 
It is satisfying that the sequence of phases along this line 
(moving right to left) is identical to that presented in 
Fig H (moving left to right) . Entirely independent theo- 
retical arguments for the sequence in Fig ^Jwere presented 
in Section [I] using methods that preserved SU(2) sym- 
metry at all stages, in contrast to those discussed in the 
present section. 

Further examination of Fig [ll] also suggests an appar- 
ently exotic circumstance under which there could be a 
second order transition between the HN and BC states. 
This happens if parameters are such that the XY+BC+I 
phase has just been shrunk to vanishing size i.e. the 
three multicritical points in Fig ^ have just coalesced 
into one, and the line of SU(2) symmetry passes through 
this coalesced point. Such a point does appear to rep- 
resent a multicritical point, requiring the tuning of more 
than one parameter, but we cannot definitively rule out 
that accessing such a point in the subspace of SU(2) sym- 
metry requires only one tuning parameter. In any event, 
this reasoning does lead to a candidate field theory for 
the critical point between the HN and BC phas es (multi- 
critical or not): it is the field theory £ in ( |5.14 ) with the 



primary fields 0i,2, b^, ^ bc and ^ i all critical, while <i> 
is massive and can be integrated out. 




FIG . 12. Schematic mean field phase diagram of £' in 
( [5.14[ ) as a function of S2 and si with si > and sbc > 0. 
The conventions for t he pha se boundaries are as in Fig§. As 
we discuss in Section V C 3 , there is no reasonable case for 
a line of SU(2) symmetry in this phase diagram, apart from 
within the XY+I and BC phases. 



3. Interplay with fractionalization 

The previous section paid particular attention to the 
role o f BC order in the phase diagram of the field theory 
£ in (5.14). This section will carry out the complemen- 
tary analysis for the case of F order. We will also ad- 
dress the issue of whether SU(2) symmetry can be found 
in such phases. We argued earlier, in Section [n| that F 
phases are not found in the class of SU(2) symmetric 
systems that have been considered so far: our discussion 
below will further support this conclusion. However, we 
remind the reader that in the following Section [v| we will 
introduce a new class of models which do have F order 
and SU(2) symmetry. 

As a prelude to our discussion, w e det ermine a section 
of the mean-field phase diagram of ( |5.14 ) as a function of 
S2 (which controls condensation of the field $ essential to 
the appearance of F order) and sj (which controls the I 
order necessary for the presence of phases with SU(2) 
symmetry). The other two 'masses', si and sbc are 
taken to be positive. The results are shown in Fig llj. 



Let us focus our attention on the F phase in Fig [12] and 
ask if it is possible for it to possess SU(2) symmetry. A 
simple argument shows that this is not possible, at least 
in the framework of the continuum models described by 
£ . From Table [V, only the field <& is condensed in the 
F phase. As was noted in Table |n|, vortices with phase 
winding +27T in the field $ carry boson charge 1/2, and 
are therefore the S z = ±1/2 spinon excitations about 
the F state. (Note that such vortices were not possible 
in the BC states of Fig [ll] and [l^ because the field 4>i.2 
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was condensed and this would be double-valued around 
a ±27r vortex in the <& field). A theory for these spinons 
follow s by setting the massive fields 0i.2 = in C in 
( 5.12 ). The resulting Lagrangian for $ is (again) s imply 



that of sca lar electrodynamics. Just as in Sections V C 1 
and V C 2 , we can dualize this to an 0(2) ip 4 field theory 



which we now express in terms of the complex scalar 
^sp. This is the field operator for the S z = ±1/2 
spinons, and is (roughly) the 'square-root' of the field 
^xy in ( |5.23| ); it is controlled by a '|* s |4 ' 



field theory 

identical in form to that for f m (5.23). Our purpose 
here was to search for SU(2) symmetry in the theory of 
\I/ S p and it is evident that this is not possible in such a 
field theory which only has a U(l) symmetry. The single 
spinon states with S z = ±1/2 are degenerate and so do 
form an SU(2) doublet; however, there is no reason, in 
the presence of the l^spl 4 interaction, for the S-matrices 
of the multi-spinon states to possess SU(2) symmetry. 

To further confirm this conclusion, let us look at states 
in the vicinity of the F state and examine whether any 
of them can have SU(2) symmetry. While the XY state 
shares a phase boundary with F in Fig [l^, note that 
the I phase does not: this prevents the XY±I state from 
being connected to the F state except across a special 
multicritical point. The general reasons for the separa- 
tion of th e I a nd F states were noted earlier at the end 
of Section 



VB 



in the discussion around (5.13) 



-this con- 
clusion also holds in the presence of BC order: one can 
show by a very similar argument that it is not possible 
for the F+BC and the F+I states to share a se cond o rder 
phase boundary. As we discussed in Section V C 1 , the 



XY±I state does posses a line of SU(2) symmetry along 
which there is HN order. We can now try to connect 
this to the F state in Fig [l2|. We clearly cannot pro- 
ceed across either the F+I state or the XY state as these 
state have anisotropic magnetic order which explicitly 
break SU(2) symmetry. The only remaining possibility 
is to access the F state directly across the multicritical 
point at si = S2 = 0. However this is also not plausi- 
ble: if the multicritical point was SU(2) symmetric, we 
would expect it to respond symmetrically to a uniaxial 
anisotropy with respect to I and XY orders, and to be 
flanked symmetrically either by I+F and XY+F states, 
or by I and XY states. This is clearly not the case in 
Fig § 

So in agreement with the arguments in Section O, we 
have found no plausible line of SU(2) symmetry in Fig [l2| 
which can extend into the F state. Consequently we have 
also not found a direct transition between the HN and F 
states. 



VI. FRACTIONALIZATION WITH SU(2) 
SYMMETRY, AND NON-COLLINEAR SPIN 
CORRELATIONS 

This section will expand upon the theory of fraction- 
alization in two-dimensional spin systems originally pre- 



sented in Refs |TJ-|T^. This formulation will exhibit the 
intimate connection between fractionalized phases and 
non-collinear spin correlations. Our treatment here will 
also incorporate insights gained from the recent analyses 
of Refs ^ ^: in particular our phase diagrams here con- 
tain a new phase (the HN+F phase) not considered in 
the original treatment. 

The treatment here will apply to frustrated S = 1/2 
antiferromagnets on the square lattice like those in fll.l| ). 
In particular, the model with first, second, and third 
neighbor exchange (the so-called J1-J2-J3 model) has 
magnetically ground states with spiral spin order, and 
fractionalized paramagnetic states could appear in its 
vicinity. Further, the same general theory should also 
apply to magnetically ordered and paramagnetic phases 
of other frustrated lattices in two dimensions, including 
the triangular lattice. The case of the triangular lat- 
tice was considered in Refs. |T7|Jl2"| : its three sublattice 
magnetically ordered state can be viewed as the limiting 
case of a spiral state on a distorted square lattice whose 
wavevector has been pinned at a special value determined 
by triangular symmetry. 

As in Refs |ljjl^, we will present the theory here using 
CP 1 degrees of freedom. In the body of the paper so 
far, we have used the 0(3) vector rij to represent the lo- 
cal orientation of the antiferromagnetic order. The close 
relationship of this previous formulation to the CP 1 ap- 
proach is discussed in Appendix |X|: a reading of this ap- 
pendix is a prerequisite to the discussion in the present 
section. The primary actor in CP 1 formulation is a com- 
plex spinor Zj a (a =j, J.) o f un it modulus (%2 a \zj a \' 2 = 1) 
which is related to by (A2), which we reproduce here 
for completeness: 
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with a = x,y, z a spin index, and a a the Pauli matrices. 
This representation has the disadvantage of introducing 
redundant degrees of freedom associated with the local 
compact U(l) gauge symmetry Zj a — > e l ^ j Zj a . However, 
as was_also the case in the treatment of Shraiman and 
Siggiacj of doped antiferromagnets, it allows for an espe- 
cially simple theory of the appearance of non-collinear, 
spiral spin correlations. Furthermore, the Zj a quanta 
carry spin S — 1/2, are therefore natural candidates for 
spinon excitations in fractionalized paramagnetic phases. 

As a prelude to writing down our model, we need to 
understand the characterization of non-collinear spin cor- 
relations in terms of the Zj a . Consider the spin correla- 
tions between two nea r ne ighbo r sit es j = 1,2. In the 
context of models like (1.4) and (A3), the optimum con- 
figuration of ni and 112 is to have them parallel to each 
other: such a state has \z* a Z2 a \ 2 = 1, and simultaneously 
£ctpz* a z 2a = 0. We would like to an energetic preference 
for the n vector to rotate by a non-zero angle in moving 
from site 1 to site 2: this can be ensured by demanding 
that \f- a (}zl a Z2a\ 2 7^ 0, while \zl a Z2a\ 2 remains non-zero. 
So, natural order parameters for non-collinear spin cor- 
relations are the two fields 
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(6.2) Suitable factors of the exponential of the compact U(l) 



where the index g will always extend over g = x,y. These 
fields are not gauge invariant and preserving gauge sym- 
metry will, of course, be essential in determining the 
structure of the effective action. Notice also that in the 
continuum limit, Q e ~ e a f3Z a d g Zf3, which is the spiral 
orderparameter used by Shraiman and Siggiau3 an d in 
Ref. [l5| The form of the order parameter chosen in ( |6.2| ) 
corresponds to spin spirals along the (1,0) and (0, 1) di- 
rections of the square lattice. Clearly, it is possible to 
have antiferromagnets in which the spiral correlations are 
oriented along the (1,±1) or other directions: there is 
a corresponding modification of the order parameter in 
these cases, but we will not present this simple general- 
ization explicitly here. 

To fully explore the phase diagram of our model, we 
find it convenient to introduce explicitly the real, Neel, 
order parameter vector Nj. This is the 0(3) symmetric 
combination of the orders ^ xy and considered in 
Section [VCl. Clearly 



gauge field Aj fl have been inserted to ensure that the 
action is invariant under the U(l) gauge transformations 



Qj e — > Qj 



J4j 
P 2i <t>j 



A, 



A, 



(6.5) 



where <f>j represents an arbitrary phase. The potentials 
in Sq and Sjy are given, as usual, by 



v i( x ) = o ( s F + u tF 2 ) 



(6.6) 



where jj = Q,N. The allowed cubic terms between the 
fields have been explicitly written out in S w . A number 
of quartic terms coupling the modulus squared of all the 
fields are also allowed but we have not written them out 
for simplicity. 

The model Zp in (6.4) is our global theory for the 



Nja ~ Zj a a ap z j0> 



(6.3) 



where a = x,y,z. Also it is evident that Nj ~ nj, but we 
prefer to use an independent symbol here as the physical 
context is different, and we do not impose a unit length 
constraint on Nj. 

We can now write down the partition function for 
quantum antiferromagnets with SU (2) symmetry. The 
expression here generalizes Z cp in (A3) of Appendix ^ 
by allowing for non-collinear spin correlations 



Zp = e / n dA ^ J n dz j a s(\ Zja \ 2 - 1) 

II d Q]Q J II CX P (-^a - S A - Sb ~ S Q - S N - S v 

je ja 

g J_y^Q. e 2iA in Q . + . + c c + V (\Q ■ | 2 ) regime is continuously connected to the strong-coupling 



phases and phase transitions of quantum antiferromag- 
nets. We believe it offers a complete theory of the univer- 
sal properties of S = 1/2 antiferromagnets in two dimen- 
sions with SU(2) symmetry on a variety of bi-partite and 
non-bi-partite lattices. Further, all its phases have the 
obvious generalization to cases where SU(2) symmetry is 
broken weakly by an easy-plane or easy-axis anisotropy: 
HN order will transform to the corresponding XY or I 
order. 

Various symmetries have played a crucial role in de- 
termining the structure of Zp, and will obviously also be 
important in the analysis of its phase diagram. The fore- 
most among thes e is the compact U(l) gauge invariance 
associated with ( |6 . 5[ ) . Along with the gauge field Aj^, 
we have two "matter" fields: the unit charge field Zj a , 
and the charge 2 field Qj e . The phases of such, a gauge 
theory were discussed by Fradkin and Shenkcia, and we 
review the essential needed results. For the pure gauge 
theory (without matter), the weak-coupling (small e 2 ) 
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S w = -WQ Q^eapZjaZj+g^e 1 ^" + c.c. 

j 



(6.4) 



As stated earlier, all sums over g extend only over x,y, 
while those over fj, extends over x, y, r. The expressions 
for the first three terms in the action were given earlier in 
(A3); these are: S z , the hopping term for the z a , Sa, the 
Maxwell term for the compact U(l) gauge field Aj^, and 

Sb, the Berry phase of the underlying 5 = 1/2 spins. 
The symmetry of the problem clearly permits a more 
anisotropic nearest neighbor coupling in Sq, but we have 
written down an isotropic form for notational simplicity. 



(large e 2 ) regime, and both of which exhibit confinement 
of unit charges. In the presence of charge 1 matter fields, 
one might expect a separate "Higgs" phase where the 
matter field has condensed: however, this charge 1 Higgs 
regime is continuously connected to the confining phase, 
and there are no intervening phase transitions (at least 
none associated with the gauge theory — other symme- 
tries could (and do) lead to phase transitions). Upon 
adding the charge 2 fields, their condensation does lead 
to a charge 2 Higgs phase, which is separated from the 
remainder of the phase diagram by a phase transition: 
unit charges are deconfined in this Higgs phase, while 
they are confined elsewhere. 

The second important symmetry is the global SU(2) 
invariance: Zj a transforms like a spinor under this sym- 
metry, while Nj transforms like a vector; all other fields 
are SU(2) singlets. This symmetry is broken in phases 
with HN or P orde r, and their properties will be discussed 
further in Section VI B below. 
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Finally, Zp in invariant under the square lattice point 
group D4. These symmetries can also be broken, leading 
to new phases. As we have seen in the previous sections 
of this paper, the Berry phase term can induce BC or- 
der, which breaks D4 symmetry In a region where Qj g 
is suppressed (and so we can neglect the presence of a 
charge 2 Higgs fields), the compact U(l) gauge theory 
is always in a confining phase, and (as we noted above) 
there can be no phase transitions associated with it; nev- 
ertheless, Berry phases can induce a transition between 
the HN and HN+BC phases (both of which are con- 
fining) associated with D4 symmetry breaking. In the 
region where Qj e is non-negligible, another type o f D4 
symmetry-breaking becomes possible. Notice from (6.2) 
that Qj e transforms like a spatial vector (the E represen- 
tation): it is odd under reflections, and its components 
map into each other under 90° rotations. Consequently, 
the Higgs phase where Qj e is condensed can also break 
D4 symmetry: we will see that the paramagnetic F phase 
can also have "bond nematic" ordenl3'C3, associated with 
the choice of condensation between Qj X 



Note 



and Qjy 

also that although Qj X — > —Qj X under a reflection along 
the y axis, condensation of Qj X does not imply loss of this 
reflection symmetry: Qj X also carries a compact U(l) 
gauge charge, and this change of sign can be undone by 
a gauge transformation. 

Careful consideration of the above symmetries, and in- 
sights from the analyses of the previous section, lead to 
sections of the proposed phase diagram for Zp in Fig [l3| 
and [l4|. We also present a summary of the character- 
izations of the various phases in terms of the fields in 
Zp (and in Z' P below) in Table fv|. We will describe vari- 
ous portions of this global phase diagram in the following 
subsections. 




FIG. 13. Section of the proposed global phase diagram of 



Zp in (6.4) as a function of sn and sq for moderately large 
values of g, qq and g^. The phases with F order are likely to 
have "bond-nematic" order over a substantial (and possibly) 
portion of their domain. Bond-nematic order corresponds to 
a preferential polarization of the spin correlations along the 
x or y axes, and spontaneously reduces the lattice symmetry 
from Di (tetragonal) to D2 (orthorhombic). The nature of 
the phase boundaries near the multi-critical points A and B 
is uncertain, and a direct second-order transition between the 
HN+F and HN states also appears possible; this is discussed 
further in Appendix [h|. Most of the P phase has coplanar 
spin order, but portions do have non-coplanar: the latter or- 
der appears in the region contiguous to the HN+F phase; the 
phase boundary between coplanar and non-coplanar spin or- 
der is not shown. This phase boundary meets that between 
the F and HN+F phases, and so the point where the P, HN 
and HN+F phases meet is tetracritical. 



A. Connection with Z2 gauge theory 



The phases of the compact U(l) gauge theory noted 
above provide the most complete description of the phase 
diagram of Zp. However, in some portions of this phase 
diagram a simpler, description may be provided by a Z2 
gauge theoryliirQEj. This reduction will also allow us 
note the connection to the work of Senthil and Fishe 
However, it is important to note that the Zi theory below 
does not provide a complete description of the phases and 
phase transitions in Figs [l3| and [l4| we will argue below 
that it does not contain the familiar Neel (HN) state, or 
its phase boundaries. 

The reduction to the Zi gauge theory is most naturally 
observed in the limit rq <C 0. Here Qj e is expected 
to be large (at least on short scales), and by the Higgs 
phenomenon, this will also suppress fluctuations in the 
flux of Ajp. However, because Qj e carries charge 2, a 
gauge flux of n is invisible to it, while the unit charge 
fields Zj a do observe such a flux. To allow only such 
fluxes, we can safely assume that e tAj ^ is pinned to ±1 
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FIG. 14. As in Fig but as a function of g and sq for sjv 
large and positive, and for moderately large values of gQ and 
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by the large values of Qj g . So we write 



(6.7) 



State 


Description 




Nj is condensed, Qj g are massive, and unit charges 
of the compact 17(1) gauge theory Zp are confined. 
As this is continuously condensed to the charge-1- 
Higgs phase, we can also view the z a as condensed. 
This phase cannot be described by the Z2 gauge 
theory Z' P . 


BC 


Nj , Qj g are massive, and unit charges of the com- 
pact U (1) gauge theory Zp (or the Z2 gauge theory 
Zp) are confined. The Berry phases induce BC or- 
der 


HN+BC 


As in the HN state, but the Berry phases induce 
BC order, and this phase also appears in Z' P . 


F 


Nj and z a are massive, while condensation of Qj e 
puts the compact U(l) gauge theory in its charge- 
2-Higgs/deconfined phase; the Z2 gauge theory is 
also deconfined and the z a quanta are liberated. 
Bond-nematic order is likely to be present over a 
substantial portion of the F state. 


F+BC 


As in the F state, but BC order is also present. 


HN+F 


As in the F state, but Nj is condensed. 


P 


Nj, Qj e are condensed, and unit charges of the 
compact U(l) gauge theory Zp (or the Z2 gauge 
theory Z' P ) are confined; consequently, as in the 
HN state, we can also view the z a as condensed. 


P+BC 


As in the P state, but BC order is also present. 



TABLE V. Summary of all the phases in Figs [TJ and 
and their physical properties. 



where Sjj+fi — ±1 is an Ising gauge field on the links 
of the d irect lattice; it is similar to the Z2 gauge field 
in (1.15), and hence to that in Ref. |8| Clearly the field 



Zj a carries unit Z 2 gauge charge, while Nj and Q 



3l> 



are 



neutral. 

Despite the quenching of the gauge charges associated 
with Qj e , we cannot entirely neglect their fluctuations. 
This is because we still have to account for their trans- 
formations under the D4 lattice symmetry and the asso- 
ciated symmetry breaking. So we introduce additional 
Ising variables, ~Kj X and 7Tj y by the correspondence 



(6.8) 



itj g also transforms under the E representation of D4. 
Ordering transitions in the TTj e will signal D4 symmetry 
breaking with the appearance of bond-nematic order. It 
is important to note that despite being Ising spin vari- 
ables, the iTj e are neutral under the Z2 gauge field Sjj + p i . 

With the above parameterizations, we can reduce Zp 
in (^J) to the following Z 2 gauge theory 

z' P = e /n^d^i 2 - 1 ) 

{s jd+il = ±l}{ir je =±l} J ja 

J Y[ dNj exp (-5; -S N - S' w - S s - 5V) 

S z = /. z ja S j:j+P- Z j+^,a + c - c - 

s 'w = - W Q £ 7r je S 3d+e £ a0 z 3aZ j+ g^ + c.c. 
3Q 

- w N ^N 3a z* a a a al3 z jl3 
3 

~ k n s ^i+/ 1 + 'f X (i - s J^+ f ) 



Ss = -1 



Sir — 1 7r je 7r j+A,e + v n Z~^, i^jx^j+p^x) ( 7r jy 7r j+A>y) 1 

(6.9) 



where S n is as defined in ( |6.4[ ) , the coupling g'q is a con- 
sequence of <7q in Zp, and K derives from the Maxwell 
term (proportional to 1/e 2 ) in (A3). We have also added 
a term proportional to , a repulsive interaction between 
the energies of the Ising variables Wj X and TTj y : this en- 
sures that when these fields order, the system will choose 
between orderings of one or the other, and the bond- 
nematic or spiral order will be polarized along the (1,0) 
or (0,1) directions. As we noted earlier, a different choice 
of variables should be made for polarization of the spiral 
order in other directions. 

The simplest application of the Z2 gauge theory, Z' P , is 
in the phases in which the SU(2) symmetry is preserved. 
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Here we may integrate out the Zj a and Nj fields and the 
action simply becomes S s + S w , along with additional 
couplings between Sjj+p, and 7T,- e which preserve gauge 
invariance and -D4 symmetry. The model S s is, of course, 
identical to th e Z 2 gauge theory obtained in the g — > 00 
limit of ( 1.15| ) and was studied in Secti on [v| it is dual 
to the fully frustrated Ising model in ( 5.10| ). It has a 
confining phase at small K with BC order, and a decon- 
fined phase at large K in which integer Z 2 charges (like 
the z a ) are liberated-this is the F phase. So the present 
model describes the F and BC phases appearing for large 
rjy in Fig [l3] and for large g in Fig |l4]. An extension of 
this model with more complicated couplings between the 
Isinggauge spins could also lead to the F+BC phase in 
Fig [14]. The ordering of the Ttj e variables leads to bond- 
nematic order: such order is already present in the BC 
phase and is therefore not significant there. On the other 
hand, in F phase the bond-nematic yields a preferential 
polarization of spin correlations along the x or y axes, 
and breaks the D4 symmetry down to D 2 . In principle, 
the bond-nematic ordering and the deconfinement transi- 
tions are independent and can occur at distinct positions. 
However, in the large N theory of Refs. ^4|JT5| these two 
transitions co-incided and the F state had bond-nematic 
order over its entire domain; we expect here that bond- 
nematic order will be present over a substantial, if not 
complete, portion of the phases with F order. 

A similar description of these two non-magnetic phases 
can also be obta ined starting from the compact U(l) 
gauge theory in (6.4). Again, we integrate out the Zj a 
and Nj and end up with the theory of a compact U(l) 
gauge field coupled to a charge 2 scalar. The latter the- 
ory was described in Ref. [l8] and possesses a transition 
between BC and F phases identical to that described by 
the Z 2 gauge theory S s . 

We now attempt to extend the above study of the Z> 
gauge theory, Z' p , to the remaining phases in Figs |l3j 
and ^ with broken SU(2) symmetry, and will quickly 
run into a problem; we will have to return to the orig- 
inal compact U(l) gauge theory, Zp, t o rectify the sit- 
uation. As we will see in Section VI B, all the magnet- 



ically ordered phases have (Nj) ^ 0. So let us assume 
that (N) = No (0,0,1). Gauge-invariant fluctuations of 
N about this direction will lead to the two required spin- 
wave modes. We now examine the gauge-dependent por- 
tions to check for the presence of other orderings. In the 
pres ence of N = (0, 0, 1)Nq and wn > 0, we see from 



(3.9) that among the CP 1 fields, the Zji fluctuations will 



be preferred from the Zji fluctuations; this corresponds to 
an easy-axis anisotropy to the CP 1 model, in contrast to 
the easy-plane anisotropy considered in Appendix |]. So 
we can safely neglect Zji , and model the other component 
by a phase variable z$\ = e l9jT . The resulting effective 
action for the angle dj\ and the Ising gauge field Sjj+p, in 



daily the phase diagram in Fig ^ to the present situation. 
Because there is background Neel order with (N) 7^ 0, 
the physical interpretation of the phases is now different. 
In particular, it is obvious that the BC phase in Fig 
corresponds to the HN+BC phase in Figs 13] and [l4|, and 
the F phase of Fig f|] maps to the HN+F phase of Fig [l^. 
The interpretation of the XY phase in Fig [7] is slightly 
more subtle: this state has a gapless Goldstone mode as- 
sociated with the XY order, and in the present situation 
this must be combined with the two spin-wave modes of 
the N order to yield a total of 3 gapless modes. This 
can only correspond to the P phase in Figs and |TJ. 
As we will see below in Section VI B 3 the N order also 



gets spatially modulated in the P phase, and this is an 
effect that the present simple analysis misses. Similarly 
the XY+BC phase found in Fig O maps here onto the 
P+BC phase of Fig |J. 

The most striking feature of the above analysis of 
phases with bro ken SU(2) symmetry of the Z 2 gauge the- 
ory, Z'p in (6.9), is that we achieved a description of all 
the phases and phase transitions in Figs 13 and [lj] except 
those associated with the HN phas e. In the above map- 
ping of its properties to Zrj(i) in (1.15), the HN phase 
requires a phase of •Z'u(i) with no broken symmetries and 
no fractionalization or topological order. No such phase 
was found in entir e an alysis of Section [v], and so we con- 
clude that Z'p in (3.E) does not possess a HN phase. So 
we are in the slightly embarrassing position of having 
a sophisticated theory for the exotic states in Figs 
and lj, but are not able to capture the simple and fa- 
miliar Neel state. This situation is rectified by retu rnin g 
to the original compact U(l ) gau ge theory Zp in (|6.4|) . 
As we will show in Section VI C| and Appendix |H], the 
compact U(l) gauge boson (unlike the Z 2 gauge field) is 
able to quench the 9-\ fluctuations without breaking any 
symmetries or inducing additional gapless modes, as is 
required in the HN state. 



B. Phases with magnetic order 

Now we describe more completely the phases in Figs [l3| 
and [li] that break the spin rotation symmetry by a 
spontaneous magnetic polarization: such states will pos- 
sess gapless spin-wave excitations. The phases in Sec- 
tions VI B 1 and VI B 3 below can also co-exist with BC 



( p.9| ) is essentially identical to the model (1.15) which was 
studied in great detail in Section [v| (the 7Tj e Ising fields 
are innocuous spectators controlling the bond-nematic 
order). So we may directly apply the earlier results, espe- 



order, and this will (obviously) not change the structure 
of the gapless modes. 



1. HN phase 

This familiar Neel state is found in the regime with 
tq > and r/v < 0. So the Q e fields are massive and can 
be neglected, while the N field is condensed. The remain- 
ing degrees of freedom, and z a , obey the standard 
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action of 2+1 dimensional compact U(l) gauge field cou- 
pled to a complex scalar of charge 1, along with an addi-. 
tional Berry phase term. This latter theory is expectedc-3 
to always be in a confining phase: the Higgs phase where 
the unit charge complex scalar z a is condensed is con- 
tinuously connected to the confining phase of the pure 
compact U(l) gauge theory. So the present phase with 
Q g massive and N condensed, also has z a condensed. We 
will present a more com plete analysis of this important 
phenomenon in Section VI C and Appendix ||: we will 
find that the Berry phases can also induce a region with 
co-existing BC order. 

This starting point of condensation of N and z a pro- 
vides a straightforward description of the fluctuations in 
the HN phase. We perform a standard Gaussian fluctu- 
ations analysis about a mean-field saddle point at which 
N and z a are non-zero, and all other fields are zero. The 
fluctuations of the gauge field are 'Higgsed' by the 
z a condensate, and are consequently massive. However, 
there are two gapless spin-wave eigenmodes which arise 
from eigenvectors made up the transverse components of 
N and z a . These are just the properties expected in the 
familiar collinear Neel state. 



2. P phase 

Let us approach the P phase from the F phase, where 
as we noted in Table the Qj g are condensed. For sim- 
plicity, let us assume that Qj X is condensed while Qj y 
is zero. By the Higgs phenomenon, the fluctuations 
are then quenched and can be neglected. From Zp we 
see then that long-wavelength action controlling the spec- 
trum of the z a spinons is then controlled by the action 

S z = J d 2 xdT[\d^z a \ 2 + A 2 \z a \ 2 

- w Q Q x e a pz a d x Z0 + c.c. ; (6.10) 

we have relaxed the length-constraint on the z a and re- 
placed it by a quadratic mass term which in duces a spin 
gap A. Note that the action for spinons in (3.1C) is ex- 
plicitly SU(2) invariant, in co ntrast to the discussion of 
the F phase in Section V C 3 . It is a simple matter to 



diagonalize the quadratic form in (3.10): we obtain two 
complex modes with dispersion 



= k 2 



A 2 ± w Q \k x Q x \) 



1/2 



(6.11) 



there are two degenerate real modes at each frequency. 
This dispersion relation is sketched in Fig p"5|a; note 
that it always has minima at non-zero wavevectors k — 
±fc m in = ±(wq\Q x \/2, 0) and the minimum energy is 

1/2 

(A 2 — Wq\Q x \ 2 /4) . For large enough Q x this mini- 
mum energy will vanish, and the z a will condense at the 
wavevectors ±fc m ; n . By a determination of the eigen- 
vectors we see that the most general condensate has the 




FIG. 15. Plot of the spinon dispersion in (a) the F phase 
(Eqn (|s7lTb), and (b) the HN+F phase (Eqn (p 7 !^). For 
larger No, the latter has minima only at k = 0. There is one 
complex (or two real) bosonic mode at each frequency. 



spatial dependence (assuming, without loss of generality, 
that Q x is real) 



id\ d\ 



-ik n 



(6.12) 



where d\^ are complex numbers, with only \d\\ 2 + |g?2| 2 
determined by higher order terms in the action. We 
should no w com pute the staggered spin polarization by 
inserting ( |6.12 ) into ( |6.3| ). This is simplified by noting 
that the first factor on the right hand side is simply pro- 
portional to a general SU(2) matrix, and will only lead 
to an overall rotation in N. So up to an arbitrary 0(3) 
rotation, we can set d\ = 1, da = 0, and then the spin- 
polarization in the condensed phase is seen to have the 
form 



N tx (sin(2fc min a;), - cos(2k min x), 0); 



(6.13) 



this is an incommensurate coplanar spiral in which N pro- 
cesses at the wavevector k = 2fc m j n , as is expected in the 
P state (recall that for the underlying spins, this wavevec- 
tor is measured relative to (tt, it)). The critical properties 
of the transition at which the z a condense (which is ex- 
pected to have 0(4) symmetry), and that of fluctuations 
in the P states can be addressed as in Ref. [12] and will 
not be discussed further here. We merely note that the 
A^ gauge fluctuations are quenched across the transi- 
tion, and so the theory can be expressed in terms of the 
z a alone; the Q are also not critical, while N is pinned 
to the z a via (KL 



3. HN+F phase 

The spectrum of spin excitat ions in the HN+F phase 
can be analyzed as in Section VI B 2 for the F phase, 
but we now have to condense N. Let us assume, as in 

JVo(Q,0,l). 



Section |VIA| , that (N) = 7V (0,0,1). Fluctuations of N 
about this state will lead, as usual, to the two familiar 
gapless spin-wave modes. For the gapped spi non s tates 
we use the same strategy as that leading to ( 3.1C| ); the 
presence of the N condensate now modifies this to 
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S z = d 2 xdr [\d, 



fl^Ct | 



■ A 2 \z a \ 2 - w Q Q x t a pz a d x zp + c.c. 



-w N N Q (\zi\* - |z;| 2 )] . 



(6.14) 



This can now be diagonalized, and we again find two 
complex modes with dispersion 



w = (k 2 + A 2 ± ^{w Q k x Q x y + (w N N y^j ' . (6 



15) 



These dispersio ns ar e sketched in Fig 15 3. Unlike (3.11), 
the expression ( 6.15| ) can either have minima at non-zero 
wavevectors or only at zero wavevector; the latter applies 
for wnN > Wq\Q x \ 2 /2. Although each frequency has 
doubly degenerate real modes, there is now no SU(2) 
symmetry in the S-matrices that describes the scattering 
of these modes. Condensation of the spinons at the non- 
zero wavevector minima of the dispersion (in the case 
where they exist) will lead again to the P sta te. Analysis 
corresponding to that in (3.12) and Q6.13 ) now shows 
that the N spin polarization is non-coplanar in the P 
phase. Further, the N condensate has components at 
wavevectors at k = ±2/c min and at k — (note again that 
for the underlying spins, these wavevectors are measured 
relative to (tt,tt)). Condensation of the spinons at the 
zero wavevector minimum (if it exists) leads to the HN 
state which has a N condensate only at k = 0. 



C. Multicritical points 

The discussion in the previous subsections, combined 
with the analyses in Section [v| has given a complete de- 
scription of the phases in Figs |l3| and |TJ, and of most 
the boundaries between them. It remains to address the 
nature of some of the multicritical points. 

The tetracritical points in Figs [l3| and 14 all appear 
to be fairly conventional. As an example, consider the 
point where the four phases, BC, HN+BC, HN+F, and 
F meet in Fig [l3|. The transitions in the vertical direction 
are described by the BC order p aram eter associated with 
the frustrated Ising model in (5.10), while that in the 
ver tical direction is controlled by the 0(3) ip 4 field theory 
in ( |5.15 ) . The two order parameters have only an quartic 
energy-energy coupling involving quadratic terms which 
are respectively invariant under D4 and 0(3) symmetries. 

It remains to sort out the physics near the multicritical 
points A and B in Fig [l3|. This appears to be a problem of 
great complexity. However, some progress can be made 
by using the observation that all four phases in the vicin- 
ity of A and B have (N) =^ 0. This permits a description 
in terms of a simplified compact U(l) gauge theory model 
which is amenable to a complete duality transformation, 
and which does not suffer from the problem of Zi gauge 
theory Z' P of omitting the HN phase. Details of the anal- 
ysis along these lines are provided in Appendix [h|. 



VII. CONCLUSIONS 

This paper has presented an extensive study of the 
paramagnetic states of two-dimensional, S = 1/2, quan- 
tum antiferromagnets with both U(l) and SU(2) spin 
rotation symmetries. We have focused on two central 
properties that characterize such states: bond-centered 
charge (BC) order {e.g. spin-Peierls order) and fraction- 
alization (F). 

Our point of departure was the 'minimal model' of 
5 = 1/2 spin systems with local antiferromagnetic spin 
correlations: the 0(3) non- linear sigma model with Berry 
phases. We were strongly influenced in this choice by 
the remarkable success of such an approach in producing 
a theory of one-dimensional antiferromagnets (reviewed 
and extended in Section IV), with all results in complete 
agreement with those obtained by bosonization methods 
and also with the solution of a variety of integrable mod- 
els. We wrote a lattice version of the minima l mo del in 
two dimensions, the partition function Z in ( |l,13|) , and 
Sections ||0 were devoted, directly or indirectly, to ad- 
dressing its properties. We argued in Section || that the 
strong-coupling limit of the SU(2)-invariant Z had BC 
order, while that of U(l)-invariant systems was shown 
in Sections |H| and [v| to also allow F states. We also 
discussed, in Sections III and [y], how the familiar mag- 
netically ordered Neel (HN) state was connected to the 
BC state in systems with SU(2) symmetry: we presented 
a number of independent arguments which supported an 
intermediate co-existence (HN+BC) state. 

We also presented an in-depth discussion of a num- 
ber of routeSj-to frar^ianalization that have appeared in 
the literatureti'ESou'EJ. These can be divided, on phys- 
ical grounds, into at least three distinct classes which we 
review below. Near the transition to magnetically or- 
dered states (or to superconducting states in class (B) 
below), these routes produce fractionalized states with 
clearly distinct low energy spectra. We do not address 
the more subtle question of whether ground states in dis- 
tinct classes can be connected to each other via some 
path in parameter space without an intervening phase 
transition — there can be considerable rearrangement be- 
tween the excited states so that one can connect states 
with very different low energy structures without a sin- 
gular change in the ground state. The routes are: 

(A) In systems with a U{\) spin symmetry, we approach 
the F state from a magnetically ordered state with XY 
order, which we view as 'supe rfluidity' in the S z = 1 spin- 
flip boson i>t — S + (sec ( 4.2l| )). This superfluidity can be 
destroyed by condensation of,-*prtices in the phase of b' , 
and the F phase is producedEwi3 by a preferential conden- 
sation of double vortices, <i>; near its transition to the XY 
phase, the F phase has a massive S z — 1/2 spinon exci- 
tation described by the complex , boso nic field operator 
^sp, as was discussed in Section VC3| . We have argued 
in Section [v| that this route to fractionalization does not 
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survive the restoration of SU(2) spin symmetry. 

(B) The second routeBi is technically similar to (A), 
but has a very different physical interpretation, and ap- 
plies to systems with SU(2) symmetry. The boson, b* is 
re- interpreted as the creation operator for Cooper pairs, 
and the XY phase is now a true superconductor. The 
condensation of double vortices (with flux hc/e) then 
leads to a fractionalized insulator with S =1/2 fermionic 
spinons (which are possibly gapless) . At half-filling, this 
route requires proximity to superconducting states, and 
so is only likely in systems with strong charge fluctua- 
tions which can allow the Mott charge gap to collapse. 

(C) The third routeEl was discussed in Section VI. Non- 



collinear spin correlations were shown to libera te bo sonic 
S = 1/2 spinons descr ibed by the actions ( 6.1C ) and 
( |6.14D , Note that ( j3TC| ) has an explicit SU(2) symme- 
try, and the number of spinon excitations is doubled with 
respect to (A). 

We now use the perspective of the above 'classification' 
to present a critical review other recent studies of S = 1/2 
antiferromagnets ip-iasip dimensions. 

The Paris groupc3~E!3 has performed large scale numer- 
ical studies on a variety of lattices. They have found ex- 
amples of BC states, in proximity to magnetically ordered 
states with collinear order, in keeping with the discussion 
in Sections [H] and III . On the analog of the pyrochlorc 
lattice in two dimensions, they have recently founctij pla- 
quette BC order, in accordance wit h t\ m, predictions of 
the sine-Gordon model of Appendix kT3pl and the quan- 
tum dimer modelEj. They also find good candidates for 
fractionalized states with a well-developed gap to all ex- 
citations: these are in regimes where the classical ground 
state is non-collinear, and so we propose that these states 
are in class (C). The case of kagome lattice remains mys- 
terious as it has a large density of very low energy singlet 
excitations. 

We|-.have already mentioned recent numerical evi- 
dence&EI for BC order in frustrated square lattice antifer- 
romagnets with first and second neighbor exchange in- 
teractions (the - / 1 -J 2 model). Capriotti and Sorella orig- 
inally claimedtil plaquette-BC order in the J1-J2 model, 
but have recently arguedE3 in favor of a state with no 
broken symmetries pear J2/J1 = 0.5. Our studies here 
and the earlier workotj do not, in principle, rule out a F 
state in the J1-J2 model, but do require a BC state to in- 
tervene (or, more exotically a HN+F state) as one moves 
towards the HN state at smaller J2. The variational ap- 
proach in Ref. |52| does not work well in the region with 
HN order; as in our theory the BC order is intimately 
connected to fluctuating HN order, we suspect that the 
approach of Ref. [5^ is also insensitive to the development 
of BC correlations. r-ir-| 

Moessner, Sondhi, and collaboratorscZrEa have-exam- 
ined generalizations of the quantum dimer modelo (dis- 
cussed in Appendix CI) to a variety of lattices. They 
present very convincing evidence for BC phases on bipar- 



tite lattices, and for an F state on the triangular lattice; 
these results are in general accord with the discussion of 
the present paper and the route (C) to fractionalization. 
The dimer model yields an appealing picture of the F 
state as Anderson's RVB phasetl, and so this agreement 
suggests that there is a natural evolution from deconfined 
phases (C) with appreciable non-collinear spin correla- 
tions to the short-range spin singlet correlations of the 
RVB state. ._. 

Kashima and ImadaEij have examined a Hubbard 
model with first and second neighbor hopping near the 
metal-insulator transition. At large enough Hubbard re- 
pulsion, U, this model reduces to the J\-Ji antiferromag- 
net, and so the corresponding region could have BC or- 
der. However the Hubbard model also induces additional 
ring-exchange terms at smaller U , and this should in- 
crease the tendency towards F states of class (C). Nearer 
the metal-insulator transition, the strong charge fluctua- 
tions could induce F states in class (BL 

The striking experimental evidenceO for a F state in 
CS2CUCI4 appears to place it very naturally in class (C), 
as ias been discussed by Chung, Marston, and McKen- 
ziec2l: the neutron scattering measurements show clear 
evidence for incommensurate, spiral spin correlations. 

It has recently become possibleES to perform large-scale 
numerical studies of a class of frustrated S = 1/2 antifer- 
romagnets with a U(l) spin symmetry. It would be very 
useful to study these in great detail, and to compare the 
phase diagrams with those in Section wL_. 

Starykh et al£3 and Moessner et alx5 have recently 
proposed a most interesting two-dimensional paramag- 
netic phase which does not fall into any of the classes 
discussed here: they have found an example of a "slid- 
ing Luttinger liquid" ("sliding icecT), in which the low 
energy excitations are well described by an infinite set of 
essentially decoupled one-dimensional antiferromagnets 
in the TL phase of Figs ^ and ||. However, it remains to 
be seen if such a phase is stable in a parameter regime 
with full SU(2) spin symmetr.y-_ 

Zaanen and collaborators^!! have introduced a Z% 
gauge theory of "sublattice parity order" in the paramag- 
netic states of a class of doped antiferromagnets. Their 
work is in a spirit similar to (C) but differs in impor- 
tant aspects: the paramagnetic state has incommensu- 
rate spin correlations, but the spins remain collinear; 
further, there do not appear to be any neutral 5=1/2 
excitations. .— . .— . 

Rantner and Wer£3 and WenEj have recently provided 
a comprehensive classification of "spin-liquid" ground 
states in two dimensions. Many of these are fractional- 
ized states which fall into our classes (B) and (C). How- 
ever, they also propose an "algebraic spin liquid" (related 
to the earlier "flux phase"Ej) in which there is a gap- 
less U(l) gauge boson strongly interacting with gapless 
fermionic spinons. We are skeptical about the stabil- 
ity of such a state. The original analysisEa dealt with a 
field theory of a non-compact U(l) gauge field, whereas 
a microscopic derivation always yields a compact U(l) 
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field. The latter possesses instantons which have been 
the center of attention in our analysis here, and which 
have been the driving forcej-behind the appearance of 
BC order. Raejtner and Wer£3 neglected the instantons, 
whereas WenEa argued for the irrelevance of the opera- 
tor creating a single instanton under suitable conditions. 
However, this irrelevance is not a sufficient criterion for 
neglecting the effects of a finite density of instantons. 
The latter effects-xan be addressed by the Kosterlitz re- 
cursion relationsEil, which examine the renormalization 
of the instanton action by integrating out a a tightly 
bound instanton/anti-instanton dipole. In spacetime di- 
mensions D > 2 such renormalizations are disruptive, 
and cannot be simply absorbed into an innocuous renor- 
malization of a "stiffness" (as is the case in D = 2 in 
the superfluid phase of the Kosterlitz-Thouless transi- 
tion) . A number of models of the statistical mechanics of 
"charges" in D = 3 have been examined in the literature 
(with interactions between opposite charges separated by 
spacetime distance R behaving like ~ — l/|i?| (Ref. 61), 
- ln(|-Rj) (Ref. ||), and - \R\ (Ref. ||)), and all results 
are the same: there is a strong renormalization induced 
by the dipoles, and this always liberates the charges at 
large scales; this liberation includes the regime where a 
naive analysis of the single vortex fugacity suggests that 
unbound charges are suppressed. The instanton gas is 
then in the Debye-Hiickel plasma screening phase, and 
there is no phase of confined charges for D > 2. It is 
this screened plasma of instantons which is responsible 
for the BC order and confined spinons we have studied 
here, and we believe it will also disrupt the algebraic spin 
liquid. 
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APPENDIX A: CP 1 FORMULATION 

We discuss here an alternative lattice model for the 
ground states of quantum antiferromagnets in t wo di - 
mensions. It is closely related to the model Z in (1.13), 
which was in turn derived from the coherent state path 
integral of quantum spin systems. The approach here 
is closer to the model of Ref. |j: it has the advantage 
of allowing a simple generalization which accounts for in- 
commensurate spin correlations and fractionalized phases 
in systems with SU(2) spin symmetry (discussed in Sec- 
tion VI), and of having a simpler large g limit. Its main 
disadvantage is that it enlarges the number of degrees of 



freedom on the lattice scale, and the connection to the 
underlying quantum spin Hamiltonian is not as direct. 

The central actors in this formulation are complex 
spinors Zj a (a =f, j) on the sites of the direct spacetime 
lattice. These obey the constraint 



Ei 



i 



on every site, and are related to the rij by 

* a 
z ja a a(3 z 3p 



i j a 



(Al) 



(A2) 



where a — x,y,z is a spin index, and the a a are the 
Pauli matrices. Trading in the rij for the Zj a introduces 
a compact U(l) gauge redundancy under which Zj a — > 
Zj a e^ j , with (j>j arbitrary. This is closely related to (|Q|), 
but will now be an explicitly realized in every term in the 
Hamiltonian. 

A companion player in the present formulation is a 
compact £7(1) gauge field A, M . This transforms under 
the above gauge transformation in the obvious manner 
Aj^ — > Aj^ — A M </>j, and the action is invariant under 
Aj^ — > Aj^ + 2ir. As we will see below, A, M is closely 
related to Ajp. 

We can now write down the new formulation of the 
partition function: 



Zcp 
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Sb = E VjAjr 
j 



(A3) 



The parallel between Z cp and Z should clear: S z replaces 

S a , Sa replaces Sa, and Sb replaces Sb- The phase 
diagram of Z cp as a function of g and e 2 should be similar 
to that of Z, as we argue below and in Appendix |]. 

It is tempting to analyze Z cp by taking the continuum 
limit of Sa valid for small e 2 : in this case Sa reduces to 
the standard continuum action for a U(l) gauge field, and 
Z cp reduces to the theory of a complex spinor, z a , cou- 
pled to a non-compact U(l) gauge field. However, this 
approach is not correct for any non-zero e 2 : it ignores 
the compactness of the U(l) gauge field, and the asso- 
ciated instantons, which confines charged fields in 2+1 
dimensions. 

A proper understanding of Z cp is achieved by an anal- 
ysis appropriate for large e 2 . A crucial identity in this 
limit is the exact expression 
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(A4) 



34 



a suitable gauge choice has to be made on the left-hand- 
side, corresponding to the choice of no in the definition of 
Aj^. Using this, we can explicitly perform the integral of 
in Z cp and map it into a model which has a structure 
similar to Z. This is most easily done at e 2 = oo, when 
the integral over Aj^ can be done independently on every 
site, without any preliminary decoupling procedure (we 
have already argued that the e 2 — > oo limit should be 
smooth, and the physics at e 2 = oo should not differ 
from that at finite e 2 ). A simple evaluation of the integral 
shows that 



Z cp {e 2 = oo) 



= Jli *i,<5(n 2 - 1) cxp (-S B - 5 r 

3 

= JJj x ^ \/2(l + n r n j+f ) j 

x n jo ( V2(1+nj - n ^ ) ) , 

H=x,y \ & / 
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where 7o,i are modified Bessel functions which are mono- 
tonically increasing functions of the their arguments, and 
increase exponentially fast at large arguments. In the 
evaluation of the integral it is useful to shift Aj^ by Aj^ 
and this leads to the Berry phase Sb expressed in terms 
of the Ajr rather than the Sb expressed in terms of the 



Aj T in (A3). For small g it is also interesting to note that 



the integral over Aju is dominated by values for which 
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although this approximation was not used in the evalu- 
ation of S n — thus the A flux is pinned to values close 
to the flux associated with n fluctuations. The explicit 
form of differs S n differs from that of S n , but the two ac- 
tions should have essentially identical physical properties 
for any finite g. Both actions are minimized when the rij 
are all parallel (a state with HN order), and both increase 
monotonically as the angle between any two neighboring 
rij is increased. The action S n does not have precisely 
the same form along the spatial and temporal links, but 
this should lead only to an innocuous and uninteresting 
renormalization of a spin-wave velocity. This establishes 
the close connection between the properties of Z and Z cp . 

Another connection can be seen in the g — > oo limit 
of Z cp . In this limit, a direct application of the dual- 
ity methods developed in Section jj establishes the exact 
result 



Z cp {g = oo) = Z h 



(A7) 



wh ere Z h is the height model partition function defined 
in ( 2.1 2j) . So all of the results for Zh discussed in Ap- 
pendixO apply directly to Z cp . 



APPENDIX B: LARGE g LIMIT WITH SU(2) 
SYMMETRY 

This Appendix will describe the evaluation o f W (a u ) 
in (2.7) for the SU(2) symmetric action S n in ( |l.9| ) for 
large g. 

One way to proceed is to accept the equivalence be- 
tween Z and Z cp just discussed, and perform the average 
not over S n , but over S z with Aj^ replaced by Aj^: this 
will lead to essentially the same results as those appear- 
ing below. 

Alternatively, we can follow a procedure that is roughly 
the inverse of that discussed towards the end of Ap- 
pendix ^| We first, formally, write down an effective 
action for Ajn fluctuations by 



-S A (Aj^) _ 



J I]dn^(n 2 - l)5{A m - A m )e' s " 



(Bl) 

The resulting action will be invariant under gauge trans- 
formation of the Aj^ and also periodic under shifts of 
Ajjj by 2ir. However, because of the complicated form of 



(1.7), this action is difficult to write down explicitly, even 
at g = oo. Nevertheless, it is evident from a "high tem- 
perature" expansion of that correlations of the flux 
decay exponentially fast i.e. 
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where E^rj) = (e^\A^A\)j, and the correlation length 
£ = 1/(2 In g). So it is reasonable to model S' A by the 
simplest lattice action which is consistent with all the 
symmetries and which reproduces ( p32j ) at long scales: 



S ' A = 2^ 12^ vxAvA ^ ~ 27T 1^) 2 



(B3) 



Now we can easily evaluate W(a fl ) by the same duality 
transformations developed in Section || and obtain 



W(ap 



exp ^-i e^; 



\aj^A u Aj\ 




exp ry ~ A K h J 

{h,} \ 3 



(B4) 



This last expression is our main result for W{a tl ) 1 and it 
clearly shows the et-Meissner nature of the large g regime; 
notice its similarity to the London-limit action for the 
electromagnetic field in a superconductor — after allow- 
ing the discrete variables a and h to be continuous, we 
see that a plays the role of the electromagnetic vector po- 
tential, an d h that of the phase of the superconducting 
order. So (B4) indicates that q-flux wil l be expelled, and 
the phase is well characterized by (2.9). More precisely 
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we can easily evaluate the expression for Z in (| 
obtain 

Z=^cxp (-/(A^j-ap) 

{%} 

where the function f(n) is given by 

00 / 2 2 \ 
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8h and 



(B5) 



(B6) 



m— — 00 



The expression (BS) for Z is very closely relat ed to that 
for the height model partition function, Zh, in (2.12); the 



latter corresponds to the choice f(n) = e n /2. More 
generally, f(n) is an even, monotonically function of n, 
and we expect that the properties of the height model 
should be quite insensitive to the detailed form of f{n). 
The expression (fB9) can be viewed as a renormalization 
of the coupling e z due to n fluctuations. 
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APPENDIX C: PROPERTIES AND MAPPINGS 
OF THE HEIGHT MODEL 

This appendi x will review the properties of the height 
model, Zh, in ( 2.12 ), and also present its equivalence 
to a number of other models of paramagnetic states of 
quantum antiferromagnets. 

It is useful to first rewrite Zh is a form closer to that 
used in earlier work. We decompose the fixed field 
into curl and divergence free parts by writing it in terms 
of new fixed fields, X* and as follows: 



4* 



(Cl) 



The values of these new fields are shown in Fig |l6| . 
Inserting (Cl) into (2.12), we can now write the height 
model in its simplest form 




where 



Hj - hj - X 3 



(C2) 



(C3) 



is the new height variable we shall work with. Notice 
that the 3^/> have dropped out, while the Xj act onl y as 
fractional offsets to the integer heights. From (C3) we 
see that the height is restricted to be an integer on of 
the four sublattices, an integer plus 1/4 on the second, 
an integer plus 1/2 on the third, and an integer plus 3/4 
on the fourth; the fractional parts of these heights are as 
shown in Fig [lqa. The steps between neighboring heights 
are always an integer plus 1/4, or an integer plus 3/4. 
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(a) 

+1/8 -1/8 +1/8 -1/8 
XXX 

-1/8 +1/8 -1/8 +1/8 
XXX 

+1/8 -1/8 +1/8 -1/8 
(b) 



FIG. 16. Specification of the non-zer o va lues of the fixed 
fields (a) Xj and (b) introduced in (Cl). The notational 
conventions are as in Fig ^. Only the \x = r components of 
yjfi are non-zero, and these are shown in (b). 
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FIG. 17. Mapping between the quantum dimer model and 
Zh . Each dimer on the direct lattice is associated with a step 
in height of 3/4 on the link of the dual lattice that crosses it. 
All other height steps are 1/4. 



1. Quantum dimer model 

The connectionil with the quantum dimer modelE£l" 
emerges upon examining the large e 2 limit of Zh- In this 
limit, the action is minimized by choosing the smallest 
possible modulus of steps between neighboring heights: 
these are either 1/4 or 3/4. It is not possible to reduce 
the number of 3/4 steps to zero, and height configura- 
tions with the minimal number of 3/4 steps have an ex- 
act one-to-one correspondence with the set of dimer cov- 
erings of the direct lattice (heights which differ merely 
by a uniform shift of all heights by the same integer are 
treated as equivalent). This correspondence is illustrated 
in Fig [l7| Temporal fluctuations of the heights between 
two such minimal action configurations correspond to 
quantum tunnelling between different dimer states, and 
hence the connection with the quantum dimer model. 



2. Coulomb plasma of instantons 

Elevate the height variable to a continuous real field by 
t he Poisson summation formula; this allows us to rewrite 
@ as 



Zh= zZj II dH ^ - y £ ( A ^%) 5 

{mj} °° 3 \ 3 

-2iriJ2mj(Hj + Xj)) 
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Now perform the Gaussian integral over the Hj and ob- 
tain the height model partition function written as a sum 
over the plasma of "charges" mj. 
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where G is the lattice Green's function in (|3.1| ). For 
large \r\, G(r) ~ 1/Mj an d hence this is a Coulomb 
plasma of charges in 2+1 dimensions. Working back- 
wards through the series of duality transformations, we 
can deduce that these charges are the 'hedgehogs' in the 
n field. Each charge also carries a Berry phase of 1, i, -1, 
—i, on t he f our dual sublattices, represented by the last 
term in (C5) — these phases are identical to the hedgehog 
Berry phases computed by HaldaneLJ by entirely different 
methods. The plasma is always in its "Debye-screened" 
phase for all values of e 2 , and so isolated hedgehog are 
"free" . As we review in the following subsection, this free 
plasma of charges always has BC order, induced by the 
Berry phases. 



3. Sine-Gordon model 

The most practical formulation of the height model is 
in terms of a frustrated sine-Gordon model. To obtain 
thi s, w e insert a phenomenological "core" action e~ EcTn 3 
in (C4), and in the limit of large E c , sum over just mj = 
0,±1 on every site. This expresses Zh as the functional 
integral over what is now the sine-Gordon field, Hj^, with 
the action 



s »° = t £ t( A ^) 2 - y c °<^( H 3 + ty)] (C6) 



The properties of the model S s g have been discussed 
at length in Ref. |5|, and so we need not reproduce the 
analysis here. The main result is that the height Hj 
is always in its "smooth" phase (and the plasma above 
is therefore in its Debye screening phase), characterized 
by a definite value of the average height H(Hj (where 
the average is over quantum fluctuations and over all 
sites), and such a phase necessarily has BC order. The 
fractional val ue of H specifies the nature of the BC order. 
As in Section V C, we introduce the BC order parameter 
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Then, from the results of Ref. |fj,|2^ we see that the 
values arg(^sc') = 7r /4, 37r/4, 5ir/4, 77r/4 correspond to 
the four BC states in Fig |l|a, the values arg^Bc) = 
0, 7r/2, it, 3tt/2 correspond to the four BC states in Fig^Jb, 
while any other value of arg(^sc') yields an eight-fold 
degenerate state involving superposition of the orders in 
Fig |l]a and b (lattice symmetries show that the states 
with $ BC — > tg C , and ty B c — ► e m7r/2 * B c (n integer) 
are all equivalent). 
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APPENDIX D: CORRELATIONS OF DEFECTS 
IN THE tp 4 FIELD THEORY 

This appendix shows how defect correlations may be 
computed in a continuum ip 4 field theory. The meth- 
ods used hpee are a generalization of an approach due 
to Halperintl We will outline the technical steps in the 
calculation, as no details were given in Ref. |3p| . 

We are interested in a field theory with action 



where D = 3 is the dimension of spacetime, and a = 
1 . . . N for a theory with O(N) symmetry. The propa- 



gator Q 1 (k) 



k at small k, while at k larger 



than some cutoff, A, Q(k) decreases sufficiently fast with 
increasing k to make the interacting theory ultraviolet 
finite. The 'mass', m§ is tuned by varying the coupling 
g in the underlying model Z . 

We consider first the case N = D = 3. Then we have 
point hedgehog defects. In the present continuum for- 
mulation, Halperin showed that the topological charge 
density can be written as 



p h (r) =5° [p„(r)]det [^Va(r)] 



(D2) 



We are interested in evaluating correlators of p/t(r) un- 
der the action S v in a perturbation theory in u. The 
main technical difficulty in such a calculation arises from 
the delta function in ( |D2| ). However this may be easily 
treated by a standard Lagrange multiplier technique. 

We illustrate the method by computing the two-point 
correlator of Ph(f) at u = 0. In the standard field the- 
oretic method, correlators of tp, and its derivatives, may 
be easily obtained by first evaluating the generating func- 
tion 



exp ( — / d D rJ(r)tp(r) 
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For simplicity of notation, we will work in this paragraph 
with a single component field (p; the gener aliza tion to 
multiple values of a is immediate. At g = 0, (D3) evalu- 
ates to 

exp Q J d D r 1 d D r 2 J(r 1 )J(r 2 )Q(r 1 - r 2 )^j , (D4) 

where Q(r) is the Fourier transform of Q{k). Now gen- 
eralize (D3) by including two delta functions of tp 
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with r 3 7^ r 4 . We write the first delta function as 
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and similarly for the secon d de lta function. We can now 
evaluate the <p averages in (D5), and then explicitly per- 
form the integrals over the two A's. Inthis manner we 
obtain our key result for the value of (DS) at g = 0: 
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where A(r) = Q 2 {r = 0) - Q 2 (r) and 
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A(r 3 - r 4 
Q(n - r 3 )Q(r 2 - ra)Q(O) 
+ Q(ri - r 4 )Q(r 2 - r 4 )Q(0) 

- Q(n - r 3 )Q(r 2 - r 4 )Q(r 3 - r 4 ) 

- Q{n - r 4 )Q{r 2 - r 3 )Q(r 3 - r 4 ) 



(D8) 



Observe that (D8) is a Gaussian in J, and so ip correla- 
tions remain Gaussian e ven in the pres ence of the delta 
function. In comparing ( [P8] ) with (D4), we notice only 
two differences. First, the overall amplitude of the cor- 
relators has been modified by the prefactor l/(27rv / A). 
Second, Q has bee n re placed by Q; notice that the ad- 
ditional terms in (D8) become subdominant as n and 
r 2 move far apart provided Q[r) is a smoothly decreas- 
ing function of r (which is the case even at the critical 
point). This latter fact indicates that at large distances 
we may as well neglect the delta functions entirely: they 
merely renormalize the correlators with an overall pref- 
actor, which is not of interest in the universal properties 
anyway. It is now a straightforward matter to compute 
the tw o-poin t co rrelator of ph by repeated applications 
of (D7) and (D8). First by sending r 3 — > n and r 4 — » r 2 
we can obtain the correlator of the gradients of ip: 



(5(<p(r))dMrM<p(0))dv<p(0)) 
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where r ^ 0, we have assumed that Q(r) is a function 
only of \r\, and the prime indicates a derivative with re- 
spect to the argu men t of Q(|r|). Now we can compute 
C/j(r), defined in ( |3.2| ), by a straightforward application 
of Wick's theorem, and we obtain Halperin's resultEl: 



C h (r^0) = 6P 1 {r)P 2 2 (r) 
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where 
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There is additional singular contribution to Ch(r) at 
r = 0, proportional to S(r), which has been computed 
by Halperin: this contribution is precisely t hat needed 
for Ch(r) to satisfy the neutrality condition (3.3). 

Our approach above can be easily generalized to de- 
velop a perturbation theory for C% (r ^ 0) in u. We 
expect that the neutrality condition (pj) will be satis- 
fied at each order in u, as it is a topological constraint 
imposed by the existence of the smooth field ip(r). 

We now have the machinery to determine the long- 
distance decay of C^(r) at the critical point of S v . The 
key observation, noted above, is that the delta func- 
tions may be neglected in determining the leading func- 
tional form of this decay (note however, that the delta- 
function related contribu tion s were essential in satisfying 
the neutrality condition (3.3), as it depended also on the 
short distance behavior of Ch(r)). The subdominance of 
the delta- function related contributions at long distances 
holds also at each order in u, and so we conclude that 
the scaling dimension of ph is given simply by 



dim[pfc(r)] = dim [det[(9 M <^ a ] 
= {9 + Vh)/2, 



(D12) 



where the second expression defines the anomalous di- 
mension rjh- The expression on the right-hand-side of 
the first equation is a local composite operator, and its 
scaling dimension can be determined in an expansion in 
e — 4 — D by standard field theoretic methods. It is 
not difficult to show that this composite operator has 
no renormalizations to two loops (although they do ap- 
pear at higher loop order); hence we immediately obtain 
(|3.6U3.5|). We also discussed in Section fl| that we do 



not expect the conservation law (3.3), which is special to 
D = 3, to disrupt the valid ity of the e expansion at e = 1. 
The scaling dimension in (D12) also specifies the singu- 
larity in the density of hedgehogs at the critical point 
of SL; for the original model Z this critical point is at 



g 2c (see Fig 0) and 
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where Pi 2 were defined in (Dll). This result corrects an 
error in Halperin's result, quoted in (6 .33) of Ref. |3(]. It 
can be checked that the result in (D15) satisfies 



d 3 rC v ^(r) = 



d»C V;llv (r) = 0. 
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The first result states that the net vorticity vanishes, 
while the second expresses the conservation of the vortex 
current. Both these properties will h old a t each order 
in u, and this allowed us to obtain ( 3.16 ). The long- 
distance behavior of correlations at the critical point 
is controlled by the scaling dimension 



dim[J M (r)] = (6 + v)/2 



dim \n >lv \d v y)\d\kp 1 \ 
4-e + ?7 + 0(e 3 ). 
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In the present situation, as was di scuss ed in Section III, 
we do expect the conservation law (D16) to play a crucial 
role in D = 3, and to impose the exact result rj v = —2 
for e = 1. 



APPENDIX E: FIELD THEORY IN ONE 
DIMENSION 

This appendix outlines the steps which relate the lat- 



tice model Z^L in fl4.18| ) to the field theory in fl4.25|) 



We first elevate the height £j to a continuous real field 
by the Poisson summation identity. The presence of the 
parity variables, ay, requires a slight generalization of this 
identity to 
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where a = 1 — 2(^(mod 2)), and / is an arbitrary func- 
tion. As is conventional, it is safe to analyze these inter- 
face models in the limit of large "core" energy E c , which 
accounts for shor t-dis tance renormalizations. Then, ap- 
plication of ( |4.25| ) to ( 418 ) shows that 



where C\(g) is a smooth function of g, finite at g = g 2c - 
The non-universal constants C± apply on the two sides 
of the transition. Notice that the singularity is very weak 
and essentially unobservable. 

We conclude this Appendix by stating the generaliza- 
tion of these results to the 0(2) case, with N = 2, D = 3; 
this applies to the spin systems with U(l) sym metry . In 
this case, the generalization of the expression— (D2) for 

the defect density is the following expressionETfor the -2e~ 4Bc cos(2ti%) - 2e~ E " tanh(JQ)sj cos(tt<£j) 
vortex current (r) : 



Z u(i) = / JJdfcexp ( £ ^ A ^ 



(E2) 



Jfi{r) = S 2 [<p a (r)] e^\dv(pi(r)d\ip 2 (r) 



(D14) 



The above methods show that at u = 0, the correlator of 
Jn{r) is given by 

C Vtltv (r) = (Mr)J„(r)) 



In terms of these variables, the BC order parameter is 
Q j,j+3: ~ £jCOs(7t</>j) and I order is determined by Sj z ~ 
A x (/>j. 

It is useful to now integr ate out some short-distance, 
high-energy fluctuations in (|E2[). We write 



fl 1 V 



(D15) 



(E3) 
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and minimize the action with respect to the massive held oo 
<j)\. This yields 



h[r s 



2tt 

ta,nh(Kd)e~ Ec sin(7r</>o(rj)) 



(E4) 



Inserting this back into (E2), and performing the gradient 
expansion for the held (f>o we get the action of the sine- 
Gordon theory in (4.25) with 



Y = 2e 



-4B C 



tanlr(A^)7r 2 e 
9 



2„-2E c 



(E5) 



The expression (E4) now shows that the I order pa- 
rameter (measured by a staggered Sj tZ ) is ~ sin(7r^>o), 
while the BC order (measured by a staggered Qj.j+s) is 

~ COs(7T0o)- 

APPENDIX F: CP 1 MODELS WITH 
EASY-PLANE ANISOTROPY 

In the body of the paper, we applied the easy-plane 



anisotrop y imp lied by (1.14) to the central partition func- 
tion Z in (1.13), and so obtained Zyny This latter model 
was studied in great detail in Section M. We also claimed 
that the physics of the phases of the SU(2) invariant Z 
was essentially identical to the CP 1 model Z cp in Ap- 
pendix [X| As a consistency check, we will apply the 
easy-plane anisotropy to Z cp here: it is then reassuring 
to hnd below properties that are very similar to those in 
Section [v|. 



Comparing (1.14) and (A2), we see that the easy- 
plane limit corresponds to \zji\ = \zji\- So we param- 
eterize Zj a = e 10 " , and then the gauge-invariant angle 
9 = 9-\ — 9 1. Inserting this in (A2), and writing the peri- 
odic couplings in Villain form, we obtain the form of Z cp 
with U(l) symmetry and easy-plane anisotropy: 



X exp Sa — Sb — S^- — Si 



j J 



1 2 

5t = 2g ( A ^' T ~ 2lTm i^ ~ Aj>) 

1 2 



(Fl) 



where 5,4 and Sb are dehned in (|A3|). We can e asily 
apply the duality methods of Sections |l| and [v] to ( |Fl| ) 
and obtain a partition function for integer valued vectors 
fields aj M ,pj M residing on the links of the dual lattice: 



J ep,U(l) 




■( a J, 



J >< 



2 ^ 

□ 



XY 






/ BC 



XY 


1 RP 




F \ 



3 g °° o g 

(a) (b) 

FIG. 18. As in Fig R but for the model Z cPiU (i) 



(F2) 



Two possible phase diagrams are shown. In (a), the BC order 
vanishes as e 2 — > 0, and the e 2 = line of the BC state is 
fractionalized. 



This model is clo sely related to (5.2) for the case of the 
on-site coupling (4.16). The connection becomes evident 
when we rewrite the last two terms in ( |F2"|) as 



' E ( e ^^vPj\ 



(F3) 



(ep^aA, 



□ 



■i/Oja) 2 



with g' = g. For g 1 — 0, the model (F2) is identical to the 
model ( |5.2| ). Adding such a 'Maxwell' term with g' = g 
does appear to be an innocuous change, and it is clear 
that the b asic physics of the class of models described by 
(|f|) and (U) should be the same. 

For the specific value g' = g, the phase di agra m of (F2) 
does have some differences from that of (5.2) in Fig0 
(however, the topological relationship of the phases re- 
mains the same) , and two possible phase diagrams are 
shown in Fig Some insight is gained by looking at 
limiting values of the parameters in (|F2|). For e 2 = oo, 
( FS ) is identical to the corresponding limit of ( ]5.2| ) : both 
models reduce to the boson-current loop model of Otterlo 
et al^B which has XY order for all values of g; this XY or- 
der corresponds to 'superfluidity' in the gauge-invariant 
angle 9. Similarly, the g — limit is also identical in the 
two models with XY order always present. A difference 
appears in the limit g = oo: now ( |F2|) becomes equiva- 
lent the height model in ( 2 . 1 2| ) whose ground state always 
has BC or der, while (5.2) reduced to the frustrated Ising 
model in (5.10) which had both F and BC ground states. 
Finally, the limit e 2 = becomes clearest in the form 
(Fl): now the gauge field A p is entirely quenched and 
we have two independent XY models associated with the 
angles These have an XY ordering transition, and 
,p-^[,l quanta become free massive spinons in the disordered 
phase: this is evidently an F phase. 

An important unanswered question is the fate of the 
, F phase found abo ve at e 2 = for non-zero values of e 2 . 
m the case of (|5.2|) we also found the F phase along the 
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line g = oo and so it was clear that F order survived The £j,j+/; frustrates the hopping of the <f>j and it is nec 



over a finite domain for e 2 > 0. This is not so clear here 



essary to properly diagonalize ( G4 ) in momentum space 



for (Fl,F2h. The compact U(l) gauge theory is confining, 



for bj^ = 0. This was performed by LFS who showed that 



and the liberation of a finite density of instantons for e > for a suitable gauge choice in (5.6), the </>j had equiva 
could immediately confine the spinons and lead to BC 
order. This possibility is reflected in Fig |lg| a. However, 
in the absence of numerical studies of we leave open 
the possibility in Fig |l8| b that the F phase survives for a 
range of value of e 2 . 

Similar considerations apply to the easy-plane CP 1 
model in one dimension, and its properties are essentially 

.(i) 



lent minima near spatial momenta (0,0) and (n, 0). We 
choose the same linear combinations of fields near these 
minima as defined by LFS. Now we are pr epar ed to take 
the continuum limit of all the terms in (Gl), and the 
result leads to (5.12). 



identical to those discussed in Section IV for Z, 



'U(l)" 



APPENDIX H: GAUGE THEORY WITH (N) / 



APPENDIX G: FIELD THEORY IN TWO 
DIMENSIONS 

This appendix will provide the missing steps between 
the model ( |5l| ) and the field theory ( |5.12| ). The method 
is a combination of that applied in-pne dimension in Ap- 
pendix |], and the analysis of LFSa. 

First, we follow the two-dimensional analog of the steps 
leading to (E2). The main change is that the real scalar 



field 4>j is replaced by a re al ve ctor field which we denote 



bj^/ir. By these methods (|5.4[ ) maps to 



Z u(1) = J JJc% M exp 



^2 E O^aA^-a) 2 



■ J2 (2e~ 4£c cos(2b jM ) + 2e- E " tanh(X d )e JJ+A oos(6j^ 



(Gl) 



where £jj+a is defined in (5.6). We can now make all the 
terms in (Gl) invariant under a non-compact U(l) sym- 
metry by introducing an angular variable i?j and mapping 



(G2) 



in (|Gl|). The representation ( |Gl|) then corresponds to a 
particular gauge choice in the non-compact U(l) gauge 
theory. The remaining analysis is then identical to that 
in LFS. First, we introduce two complex scalars, the field 
operators for single vortices, </>j ~ e t§j , and the field op- 
erator for double vortices <I>j ~ e 2l ^ j . Then, the second 



term in (Gl) is proportional to a simple hopping term 



for <% 



E 



(G3) 



0, 



This $j quanta clearly have a minimum near k 
and so a naive continuum limit may be taken on (G3). 
In contr ast, the dispersion of </>j is provided by the last 
term in (Gl), and is 



E 



(G4) 



This appendix will study the vicinity of the points A 
and B in Fig [ll| We wi ll use the same strategy as that 
followed in Section VI A but wil l apply it to the compact 
U(l) gauge theory Zp in (6.4) rather than to the Zi 
gauge theory Z' P in fl6.9|). 

T he init ial strategy is the same as that followed in Sec- 
tion VIA for the Z^ gauge theory. We quench the field 
Nj at its condensate Ao(0, 0, 1) and neglect its sp in- wave 
fluctuations (as in the discussion in Section III A , the lat- 
ter will lead to weak additional power-law tails to the 
interactions considered below and do not alter the basic 
physical properties). This allows us to preferentially fo- 
cus on Zjf and neglect Zj\\ we write the former in terms 
of an angular variable as zjf = e l6 n . Note that this is an 
easy-axis anisotropy to the CP 1 fields, in contrast to the 
easy-plane limit considered in Appendix [f]. For the Qj B 
degrees of freedom in (6.4), we assume, for simplicity, 
that the spiral correlations have preferentially polarized 



in the x direction: so we write Q, 



„i8j 



and neglect 



Qjy With this transformations and (innocuous) approx- 
imations we can obtain from Zp the following effective 
action for the angular variables 9ji, 9jq and the com- 
pact U(l) gauge field Aj^ (all terms have been written 
in a Villain form to facilitate duality transforms): 



E 



jQ 



x exp Sa — Sb — Si — S, 

1 2 



(HI) 



where Sa and Sb are defined in (|A3|), and S| is defined 
in (Fl). Note that this action differs in structure from 



( Fjj ) only in the prefactor of 2 for the gauge field in Sq; 
this double charge is, of course, crucial to the structure 
explored here. The duality methods de veloped in Sec- 
tion H and are easily applicable to (HI): we obtain the 
following action for integer valued vectors fields aj^pj^ 
residing on the links of the dual lattice 

z N = ex p(-y( a J>^ a jM) 2_ ^E( £ ^ aA ^ a ) 2 
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FIG. 19. Phase diagram of Zn in (HI) and (H2) as a func- 
tion of the three couplings g, gg, and e 2 . The labels of the 
states account for the background N condensate needed in 
the derivation of Zp, and are appropriate to the phases of Zp 
in ( |6.4| ) . The figure is to be visualized by imagining that these 
couplings run along the axes of a three-dimensional cube; the 
sides of the cube have been opened out flat into the plane of 
the paper. The arrows indicate that the associated coupling 
runs from to oo along its direction. 



2 ^ 

□ 



(H2) 



Some physical properties of Zn will be mo re tr ansparent 
in this dual formulation. Note also that (H2) has only 
positive weights and so should be amenable to Monte 
Carlo simulations. 

We are interested here in the phase diagram of Zp 
as a function of the three couplings g, gQ, and e 2 . Much 
insight can be gained by looking at various limiting cases: 
the results of such an analysis of the phase diagram are 
displayed in Fig ^[ Our primary is interest is in the 
phases at center of the cube shown in Fig [l^ where the 
states found at the boundaries will meet and compete. A 
complete understanding of this strongly coupled r egio n 
will probably require Monte Carlo simulations of (H2), 
which we have not undertaken. 

The identification of the phases in Fig [ll] accounts for 
the presence of the background No condensate, a nd s o 
the phase labels are appropriate to the theory Zp in (|6.4[). 
However, taken by itself, Zn is a model of a compact 
U(l) gauge theory coupled to two angular variables, and 



its phases can also be described in the language of the 
models of Section 0. Z N possesses a gauge-invariant an- 
gular variable, 0j — Oq/2, w hich c an acquire "superfluid" 
or XY order; as in Section VIA this XY order maps to 
the P phase a fter a ccounting for the No condensate. Also 
as in Section VI A , the gauge degrees of freedom can in- 
duce BC or F phases, whi ch b ecome HN+BC and HN+F 
phases in the context of ( <5A ) and Fig |ll| 

The relationship between Zn and the models of Sec- 
tion ^ can be made explicit by examining the limit 
gq — > 0, where the charge 2 scalar e t0Q condenses. Then 
it is evident from ( |H2| ) that Zn at gQ = is exactly 
equivalent to ( 5.2Q w hich is in turn related by duality to 



the Zi gauge theory ( 1.15 ). The results of Section! 
lead to the phase diag ram a t the bottom of Fig HE . 
had argued in Section VI A that Z' P mapped onto -Zu(i) 
in the region with (N) ^ 0; as Z' P was also obtained from 
Zp in the limit of a strong charge 2 condensate, the con- 
nection between Zn and the Z 2 gauge theory (1.15) for 
small gQ is just as expected. 



We also discussed in Section |VI A that Z' p , and hence 



the small gQ limit of Zn, misses a crucial piece of physics: 
it is not able to describe the HN phase. The present 
model Zn is able to repair this key deficiency in the op- 
posite limit of large gQ . The most important property of 
Zn in this limit is that it possesses a phase without any 
broken symmetries or fractionalization, and this leads to 
the HN phase in Fig [l9]. In this respect the compact U(l) 
gauge the ory Z n is far removed from the Z2 gauge theory 



-U(l) 



in ( 1.15 ) which did not have any such featureless 
phase; this phase appears in the region corresponding to 
that with XY order in the Z2 gauge theory-unlike the dis- 
crete Z% gauge fields, the 17(1) are able to quench the su- 
perfluid order by "Higgsing" it into a massive phase. The 
appearance of HN order is, of course, satisfactory, and 
was the reason we stated our preference for Zp over Z' P . 
To demonstrate the existence of this featureless phase, 
consider the limit of large gQ in (H2): there we can set 



e^xA^pjx — 0, and Zn reduces to precisely the model 
studied in Ref. || by Monte Carlo simulations-this is a 
model of a compact U(l) gauge field coupled to a 
single charge-1 scalar. In its simplest form, such a model 
has only a single featureless phase continuously connect- 
ing the confining and Higgs regions; however, with the 
Berry phases present here a BC ordering transition can 
occur. The resulting phase diagram is shown at the top 
of Fig H 

Moving on, we can also look at the limit of large g 
exhibited in Fig [l9]. Here, the charge-1 scalar can be 
ignored, and we have the theory of a charge-2 scalar cou- 
pled to a compact U(l) gauge field with Berry phases. 
This was examined in Ref. |l8| and the resulting phase 
diagram appears on the right of Fig Its phase tran- 
sition involving competition between BC and F order is 
described by the fully frustrated Ising model in (5.10), 
and this last mapping becomes exact for small gQ. 

Finally, we consider the limit of small g appearing on 
the left of Fig [19. Now the charge-1 scalar completely 
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quenches the gauge field, and the theory becomes that of 
an ordinary three-dimensional XY model; the XY order- 
ing transition is that between the P and HN phases. 
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